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1. INTRODUCTION 


Previous studies established categories for primary and auxiliary 
electric propulsion in near-earth missions (References 1 and 2, respec- 
tively). In addition. Reference 3 characterized missions to be deployed 
from the Space Transportation System (STS) in low earth orbit. This study 
extended previous results to encompass integrated electric propulsion sys- 
tems. The objectives of this effort were to establish the nature and 
characteristics of electric propulsion systems (1) that can provide the 
propulsion required for all phases of the LEO to GEO mission, starting 
with the spacecraft assembled in low earth orbit, and continuing through 
spacecraft disposal after conclusion of on-orbit operations, and (2) whose 
characteristics have minimum sensitivity to changes in mission require- 
ments, thereby assuring as wide an applicability of the systems as 
possible. 

Curing the study effort, a low thrust computer program was developed 
and implemented for parametric analysis of electric propulsion systems. 

The gravitational equation was solved for determination of thrusting time 
for orbit boosting as a function of payload mass ratio. Th electric pro- 
pulsion system mass included thrusters, power source, power conditioning, 
gimbals, propellant, propellant storage and distribution, striicture and 
thermal control for propulsion equipment. The payload included all other 
mass deployed from the Shuttle Orbiter, In this manner, the power source 
for propulsion, which was sized for orbit transfer, was available for pay- 
load power sharing on-orbit. 

The principal parameters of interest in the analysis included thruster 
efficiency, specific impulse, specific mass, specific power, specific 
area, power fraction, and propellant fraction in addition to mass ratios 
and maneuver times. For propulsive functions other than orbit boosting, 
the velocity increment (av) characterizing the maneuver was employed, and 
the rocket equation was implemented for analysis. The propulsion func- 
tions studied included orbit boosting, inclination change, attitude con- 
trol, stationkeeping, repositioning, disposal, and power sharing on-orbit. 

In order to compute transfer times for tradeoff comparisons, it was 
necessary to add the ef 'ects of atmospheric drag, occultation (shadowed 
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portions of each orbit), and »ad’ation degradation to the climb equations. 
After these effects had been taken into account, it was determined that 
electric propulsion could be initiated at Shuttle's maximum cargo carry- 
ing altitude. It was further established that an optimum solar cell cover 
glass thickness could be used for orbit transfer. Strategies were devel- 
oped for orienting solar arrays along the velocity vector to provide mini- 
mum drag area at low altitudes, and for rolling the spacecraft to achieve 
optimum array illumination during ascent. Plane change was initiated at 
higher altitudes to minimize radiation degradation. 

The general mission parameter set investigated included payload frac- 
tions ranging from 20 to 50v, with 10 years operational lifetime at GEO. 
TransfeM' time calculations were performed for thruster ef f icier cies from 
20 to 100'!. and specific impulse up to 5000 seconds. From the resulting 
tradeoff curves, and comparison with existing technology, recommendations 
for future technology efforts were made. 

Integr-ated propulsion analysis was also performed to identify advan- 
tages of using electric propulsion for functions other than orbit trans- 
fer alone. These functions can be accomplished by using the electric 
propulsion system as an integrated part of the space vehicle rather than 
as a separate orbit transfer vehicle only. A generic spacecraft config- 
uration, that is compatible with a single Shuttle launch, was chosen for 
analysis. Thruster arrangements, thrusting modes, and functional alloca- 
tions were studied for all mission phases. Solar pressure effects on GEO 
stationkeeping requirements were found to be significant compared to 
gravitational requirements. Comparisons were made with chemical auxiliary 
propulsion systems to quantify the advantages of integrated propulsion in 
terms of weight savings and concomittant transportation cost savings. 

During the integrated propulsion analysis, it was determined that 
gravity gradient torques could be used during ascent phases of the mis- 
sion to assist in rolling the spacecraft to orient the solar array for 
optimum illumination. This was the general strategy employed during the 
s tudy . 

A summary of the technical effort is contained in Volume I of the 
final report. Details of the propulsion analysis are presented in this 
volume. 
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ORBIT TRANSFER ANALYSIS 


o 

C « 


2.1 CALCULATION APPROACH 

"P 'imary propulsion" often refers to payload deployment only, with 
"secondary propulsion" covering everything else. There are two shortcom- 
ings of this viewpoint. One is that deployment is generally thought of 
as the sequence of events that starts with lift-off and ends with orbit 
insertion, with all propulsion events being done impulsively. This rele- 
gates low-thrust propulsion, i.e., continuous as opposed to impulsive, to 
orbit and attitude maintenance; and this may not be so. For large and/or 
flexible spacecraft, some of the deployment propulsion may necessarily be 
low thrust, whether chemical or electric. The other shortcoming is that 
it places orbit changes, both plane and altitude changes, in the same 
category as stationkeeping and attitude control. Thus a misleading dis- 
tinction is drawn between a "secondary propulsion" maneuver and a "primary 
propulsion" one in terms of intent rather than the propulsion required 
for the maneuver. 

A simple way around this confusion is to regard all orbital maneuvers 
as prime. y propulsion. Then propulsion used in deployment, disposal, 
return, retrieval and maneuvers that involve altitude and/or plane changes 
that may be needed for some missions, are all included in "primary pro- 
pulsion"; while "secondary propulsion" encompasses attitude control, 
stationkeeping and docking requirements. 

From ‘‘is point of view, the orbit transfer portion of the present 
study ; an example of low thrust primary propulsion that could be done 
either chemically or electrically. The other low thrust primary propul- 
ion functions mentioned above could also be done with either chemical or 
electric propulsion. 

2.1.1 Cl imb tq^uations 

T’U' calculation of transfer time is done by integrating the radial 
. *te equation that describes the effect of essentially continuous tangen- 
tial thrust applied to a quasi-circular orbit; i.e., it does not involve 
the "characteristic velocity" av. The relations are derived below in 
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parametric form for both electric and chemical low thrust to clarify the 
similarities and differences of the propulsion types. 

From equilibrium of radial forces. 



(gravitational equation) 


so that 


V 


2 



( 1 ) 


where v and r are orbit velocity and radius and m and M are spacecraft and 
Earth mass. G is the universal gravitation constant; g, used later, is 
the acceleration of gravity. Using the Earth radius, r in meters is the 
altitude H plus the Earth radius 6.37 x 10®; in MKS units, MG * 4.00 


The energy relations are: 


Kinetic, 



potential , 


and total , 



E « U + K 


Qi^n 

IF 


(3) 

(4) 


The net tangential force, thrust minus drag, increases the total energy 
according to the relation 


(T - D) V 


dt 


(5) 


Thus, from Equations (1), (4) and (5) 


T-D » 



4 



or 



Equation (6) applies to either chemical or electric propulsion. The 
subsequent developments are different for these cases but are done in a 
parallel way for ease of interpretation. The essential difference arises 
from the fact that dry propulsion mass for electric propulsion M^p is more 
conveniently related to power than to thrust, while the reverse is true 
for chemical propulsion. 

2 . 1 . 1 . 1 Low Thrust Electric 


Thrust is given by. 



(7) 

where P is power delivered to the thruster and n and I^p 
efficiency and specific impulse. 

are thruster 

Drag is 


D=^(^)v^P 

(8) 

where A is the total drag area including the payload Ap. 

and the electric 

propulsion system Arp, v is flight velocity, p is atmospheric density and 

Cjj is the drag coefficient related to the area A. Thus 


A _ A^p ApL Mp,^ Mgp 

P ■ P MpL M^p P 

(9) 

With the definitions 


Mrp 

Or = p — {power related specific mass) 
^ 0 

(10) 

3r * (specific power) 

t «^p 

(11) 

p 

® ^ (power fraction) 

f\ 

(12) 


0 
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Equation (9) becomes 


ORlG.i'''*- ^ 

OF POOR • 


, A _ 1 , E 

V!!ep 

Vl *'pl 

where M^p is electric propulsion dry mass and Mpj^ is payload mass. 
Relating the mass flow rate to power gives, 

2nP 


( 13 ) 


m = - 


(gijp) 


I 


(14) 


and the spacecraft mass m can be expressed in terms of payload mass in the 
form. 


m = M 


PL 


1 + 


'EP 

PL 


M 



(15) 


where {Mpp)^ is the initial propellant mass and Xpp is propellant fraction. 
By combining Equations (10) and (12), 

p _ ^P^EP 

“E 

and using Equations (1), (7), (8), (13), (14) and (15) with Equation (6), 
the following climb equation is obtained: 



The propellant fraction comes from the relation 


Xpp - 
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OKlGir'iMl. FA'i'. '■'j 
OF POOR QUALITY 

By using Equations (10), (2) and (14) it is expressed in terms of the power 
fraction integral. 


^PP ~ ^ ' 


_ilL 




f 

Jo 

rPL Jo 


EP 

PL 


Xpdt 


(17) 


2. 1.1, 2 Low Thrust Chemical 


A parallel development is now done for the chemical case. Instead of 
relating m and thrust to power, they are related directly to each othe: , 


- m = 


gi 


SP 


(18) 


Drag is also related to thrust rather than power. 



and the specific area, in terms of thrust, is 


A _ \h 
T ■ T 



^PL 

^TH ^ 


( 20 ) 


where the subscript TH has replaced EP to specify the dry propulsion sys- 
tem area and mass. Also, for the chemical case, 






(thrust related specific mass) 
(specific thrust) 


so that Equation (20) becomes 


A _ 1_ 
T ■ Br 


Ot/- 


^PL 

ApL 


(21) 

( 22 ) 


(23) 
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m 


' i 



I 

OR'C-'l' - - 

OF POOR Q'J*V«: ; 


Using Equations (1), (18), (19), (21), (22), and (23), Equation (6) 
becomes 



For this case, the propellant fraction involves thrust related param 
eters. From Equations (18) and (21) 

m = . T ^ _ ^fH ^ _ _%i_ 

j 55th “c 

y‘SP T 

The propellant fraction is then 



2.1.2 aV Maneuvers 

For mission requirements expressed as a velocity change aV. the 
rocket equation is used to find the propulsion mass ratio and thrust time. 

2.1.2. Low Thrust Electric 

Using Equation (18) for thrust, the acceleration relation can be 
written , 

■ 9^SP “ 3t 

After integration between limits corresponding to total propellant expulsion, 
this becomes, 

Av 

l!lo „^^SP 

m ■ ® (rocket equation) 

8 
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or, by using Equation (15), 


^PL ^ ^EP ^ ^PPo 
MpL ^ M^p 

The propulsion mass ratio is then. 


_AV__ 
= e 9^SP 


S -(^jI ■ • ft • ■) 


The thrust time is 


( 26 ) 


(27) 


t = - 


'PPo 


(28) 


jith Equations (10) and (14) this becomes 


(glcp) 


t = - 


SP' "‘E 


2n Af 


l^p) 

\hp/o 


From Equation (26) 


\ *^EP ^ Wp 


_i‘V_ 

+ nie ^hp - 11 


and the thrust time becomes 


' 2^p 

2. 1.2.2 Low Thrust - Chemical 


/^k 

\«EP 


,i V_ 

+ llle ^^SP - 1 


(29) 


(30) 


(31) 


Equation (27) remains the same for this case, except for the jse of 
the subscript TH to designate the dry thruster mass. With this change. 
Equation (27) becomes 


MpL 


VVlA 


. -y._ 

= e ^^P 


ft •')■ 


(32) 
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Equation (28) also remains the same as for the electric thruster, but 
Equation (18) is used for the mass flow rate and Equation (21) is "ced for 
oj.. Thus in place of Equation (29), the thrust time becomes 

Combining this with Equation (30) gives 

t iVe - 

2.1.3 Comparison of Performance Relations 

The low thrust equations developed above are compared in Table 1 so 
that their similarities and differences can be more easily seen. 


Or 


SP 


V”th 



Except for the power fraction Xp, and the different specific mass 
parameters and a^., the denominators of the orbit raising equations are 
the same. Also the drag terms in the numerator are essentially the same 
and become unimportant as p decreases rapidly early in the climb. The 
remaining two terms differ by the factor 2 n/gl 5 p in the electric equation. 
Thus, for equal electric propulsion mass ratios (dry thruster plus initial 
propellant mass divided by payload mass) the climb rates have the follow- 
ing relation: 



(35) 


Also a comparison of time for the same aV maneuver, shows that 

glcp 

(t)^ = “2^ it)^ (36) 

which has the same meaning as Equation (35). For il lustration, the com- 
parison can be made more specific by using typical values for the 
parameters. 

The value s 0.024 kg/U has been used for the electric thruster 
with an uncertainty range up or down of about a factor of 2. Similarly 
qq « 0.34 kg/N is a reasonable estimate for the chemical thruster, with 


ik 


10 



Table 1. Low Thrjst Chemical and Electric Propulsion Equation SuiUTwry 

ELECTRIC CHEMICAL 
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about the same uncertainty range. Using these values in Equation (36), 
with g = 9.8 and assuming 2nXp = 1.0 

(t)^ = 0.7 I5P (t)j, (37) 

Equation (37) shows that, if a comparisor. were made on the basis of 
keeping the same propul sion-to-payload mass racio, and the same I^p, 
electric propulsion would take longer than chemical propulsion by the 
iactor of 0.7 I^p. If the I^p used for the comparison were 200 seconds, 
it would take 140 times as long for .riectric propulsion to do the same 
job; and for 4000 seconds I^p it would take 2800 times as long as for 
chemical propulsion. This is the other side of the coin that must be 
kept in view when comparing the payload advantage of tiectric propulsion. 

The payload advantage can be seen by comparing the mass ratio equa- 
tion in Table 1. The conclusion is simply that the mass ratios are the 
same for both propulsion modes when the dry thruster to payload ratios, 
and the Ijp’s, are the same. The payload advantage of electric propulsion 
derives from its higher I^p capability. 

This effect is shown graphically in Figure 1 In the lower set of 
graphs the ratio of total propiilsicn mass to payload mass is plotted 
against aV for five different values of the thruster-to-paylcad mass ratio. 
All of the individual mass ratios can be found by using these plots. 

The improvement in payload-to-propulsion mass ratio is apparent, but 
it should be noctd tnat the effect is not linear with I^p. It is a 
familiar result that propellant required for a given impulse is inversely 
proportion.'! to I^p. This intuitive idea is sometimes extended to indi- 
cate that payload mass for a given aV varies with I^p in the same way. 
Figure 1 not only shows this is not true, but also that the Ijp advantage 
incre<’ises markedly with aV and by a different ratio depending on the 
thrus :er mass ratio. In general, the AV needed for a significant Ijp 
advantage to appear increases as the thruster mass ratio is Increased. 

For low thruster mass ratios, the payload ratio increases faster than 
the Ijp ratio (for large aV) while it is considerably less than the I^p 
ratio for high thruster mass ratios. For example, at aV ■ 3000 iVsec, 
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the payload ratio is 1.4 for a 10-fold I^p increase (from 200 to 2000 
seconds) at a thruster mass ratio of 0.1, while it is only 6 for the same 
conditions at a thruster mass ratio of 1.0. 

Because the trruster mass ratio is constant in each of these plots, 
the rise in propulsion mass ratio with increasing aV comes from increased 
propellant mass. The "base-line" in the plots corresponds to the thruster 
mass ratio. The rising base line increases the propulsion mass ratio 
relative to the payload mass, but it also reduces the trip time as shown 
in the center row of plots in Figure 1 for electric propulsion and in the 
upper row for chemical propulsion. 

For electric oropulsion, the trip time is shown as the product nt, as 
indicated in Table 1. Consequently the actual time varies as l/n- This 
has a substantial effect. For example a 40% increase (e.g., n = 40% to 
56%) riisults in a 30% reduction in trip time. 

The difference between the linear and quadratic I^p terms in Table 1 
causes a qualitatively different effect of on trip time for the chemi- 
cal and electric thrusters. For the former, time monotonically decreases 
with increasing I^p while for the latter it generally increases. However 
at high aV there is an optimum I^p (for minimum trip time). At aV = 10^ 
m/sec, this is about 800 seconds for ail values of thruster mass ratio. 

At aV = 7 X 10 m/sec, the optimum I^p has decreased to about 500 seconds. 

The description above illustrates some of the important features of 
chemical and electric low thrust prooulsion. More generally, it is a 
graphical representation of the Table 1 equations for chemical and elec- 
tric propulsion maneuvers that can be prescribed in terms of AV require- 
ment. 

2.2 SPACECRAFT PARAMETERS 

The spacecraft parameters for the low-thrust, electric, orbit trans- 
fer vehicle appear in Equation (16). They are discussed individually here. 

2.2.1 Thruster Efficiency 

As mentioned earlier, time varies inversely with thruster efficiency. 
This is shown explicitly in Table 1 for aV maneuvers and it is also true 
for orbit transfers in general. For this reason, considerable technology 
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development effort is devoted to thruster efficiency improvement. For the 
purpose of this study, the three efficiency functions shown in Figure 2 
are defined. The Case I and Case II efficiency functions are lower and 
upper bounds respectively of present-day high performance electric 
thrusters. The Case III function is an estimated upper limit that may be 
approached with further development. 


2.2.2 Specific Mass 

The specific mass parameter is defined by Equation (10). It is the 
ratio of the total power-related mass to the beginning-of-life solar array 
power at normal sun incidence. The power-related masses comprise the total 
solar array mass (blanket, boom, supports, etc.) and the mass of the power 
conditioning and thruster hardware including the empty propellant tank 
and thermal control system. To put it more simply, it is the total 
electric propulsion system mass minus the propellant mass. 

Table 2 lists estimated component masses for two types of advanced 
argon electric propulsion systems. One is the pulsed inductive plasma 
thruster and the other is the argon ion thruster. 

The power processor is clearly the most important contributor. It is 
followed by the considerably smaller value of the thermal control system. 
The sum of these two components constitutes about 90% of the total elec- 
tric propulsion system mass (excluding power source). 

The solar power specific mass is affected by the cover thickness. A 
thin, advanced technology, silicon solar cell having a cell thickness of 
2 mils and cover thicknesses equivalent co 1 mil of silica has been esti- 
mated to provide an array specific power of - 0.005 kg/watt (Refer- 
ence 4). A linear variation of with total thickness of cell plus covers 
was assumed for cover thicknesses of 1, 3 and 6 mils of silica on each 
side of the cell. These values are plotted in Figure 3. The SEPS solar 
cell, also shown in this figure, was used oS ? basis for extrapolating 
the solar cell array specific mass to a total thickness of 26 mils corre- 
sponding to the 2-mil cell with 12-mil covers on each side. 


The above estimates of electric propulsion specific mass are summar- 
ized in Table 3 where the total electric propulsion specific mass is the 
sum of the thruster and solar array values: + o^. 
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Table 2. Estimated Specific Mass for Advanced 
Argon Thruster Systems 


Component 

Plasma 

(kg/kW) 

Ion 

(kg/kW) 

Power Processor 

4.60 

13.61 

Thruster 

0.36 

0.18 

Gimbals and Support 

0.07 

0.03 

Thermal Control System 

1.54 

2.53 

Propellant System (Dry) 

0.37 

0.92 

Total 

6.94 

17.27 


Table 3. Electric Propulsion Specific Mass Estimates 


Thruster Type 

“TH 

(k9/W) 

Cell Cover 
Thickness 
(mils) 

“s 

(kg/W) 

“E 

(kg/W) 

Plasma 

0.007 

1 

0.0050 

0.012 



3 

0.0059 

0.013 



6 

0.0085 

0.016 



12 

0.0250 

0.032 

Ion 

0.017 

1 

0.0050 

0.022 



3 

0.0059 

0.023 



6 

0.0085 

0.026 



12 

0.0250 

0.042 


2.2.3 Specific Power 

The specific power is defined by Fouation (11). It is the ratio of 
the begiiining-of-1 ife power to the solar cell array frontal area. The 
projected value of ISO (w/m ) from Reference 4 was used for this param- 
eter. The local power per unit mass at any instant during the climb is 
the quotient of the local power fraction Xp with the specific mass 
The power fraction is the product of the individual power fractions for 
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occultation, reduced array illumination that can occur for some spacecraft 
orientations and solar cell degradation due to radiation. 

2.2.4 Specific Area 

The specific area ^pj^/Ap^ appears in the second term of the bracketed 
expression in the numerator of Equation (16). It is the ratio of the pay- 

load mass to the payload drag area. For the orbit transfer calculation a 

2 2 
nominal value of 100 kg/m was used. A lower value, 30 kg/m , was used 

for low altitude flight calculations. 

2.2.5 Propulsion Mass Ratio 

The propulsion mass ratio M^p/Mpj^ appears in Equations (15) and (16). 
It is the ratio of the electric propulsion system dry mass to the payload 
mass. This is the input parameter for Equation (16) that determines the 
rate of climb. 

2.2.6 Propellant Mass Ratio 

This ratio, (Mpp/Mp|^), also appears in Equations (15) and (16). It 
is the ratio of the propellant mass used for orbit transfer to the payload 
mass. This ratio determines the end-of -climb altitude. The propellant 
ratio at any instant during the climb is the product of the initial value 
of this ratio with the propellant fraction calculated by continuous inte- 
gration of Equation (17). The climb is finished when this function 
reaches zero. 

2.2.7 Auxiliary Propulsion Mass Ratio 

The auxiliary propulsion mass is also expressed dimensionlessly in 
terms of payload mass, AMpp/Mpj^. This ratio does not appear in the climb 
equations derived earlier. They describe orbit transfer for the special 
case of no auxiliary propulsion requirement. The procedure described here 
is used for the general case which does require auxiliary propulsion. 

The following assumptions are made for this analysis: 

• The auxiliary propulsion dry mass is ignorable compared to the 

main thruster system dry mass (of the order of 1 %). 

9 Attitude control functions during orbit transfer (including 
inertia wheel unloading) are done by gimballing the main 

thrusters, at effectively no cost in additional propellant mass. 
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• On-orbit attitude control propellant mass is ignorable compared 
to stationkeeping propellant mass. 

• Attitude control during eclipse is done by an inertia wheel. 

With these assumptions, the integrated electric propulsion system 
masses are those due to the main thrusters M^p, the propellant used in 
orbit transfer Mpp, and the propellant used for on-orbit stationkeeping 
and disposal AMpp. 

The total electric propulsion mass ratio is therefore: 

/ **EPs\ (35 

V"PL y 

where the subscript EPS refers to the electric propulsion system (dry 
thruster system including solar power and power conditioning) plus propel- 
lant; and the subscript T indicates the total system mass ratio when aux- 
iliary propulsion is used. With the above assumptions, AMpp is the aux- 
iliary propulsion mass. 

The orbit-raising calculation is done as described above except that, 
in this case, it is necessary to distinguish between the actual payload 
mass which is used as the reference mass in Equation (35) and the refer- 
ence mass Mpi^ + AMpp that is used in the computer program for the orbit 
transfer calculation. This is done with the following notation for 
"electric mass ratio" 


MRE 


M 


EP 


MpL + AM, 


PP 


and "propellant mass ratio" 


MRP 


M, 


PP 


MpL + AM, 


PP 


(36A) 


(36B) 


The rocket equation is used to express AMpp in terms of the character- 
istic velocity increment aV associated with the total on-station and dis- 
posal impulse requirements. Thus 


^PL ^EP ^^PP 

MpL + M^p 


Exp 


AV 

gisp 


(37) 
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The state of the art auxiliary propulsion I^p = 3000 seconds is used for 
all calculations. 


By dividing the numerator and denominator by Mp|^, this becomes, after 
rearrangement. 



From Equation (36A) 



so that Equation (38) becomes 



Similarly to Equation (36A), Equation (36B) can be rearranged to give 




M 


MRP 


With the use of Equations (39) and (41), Equation (35) becomes 


("epsV 


MRE + MRP + (MRE + MRP 


• "(a 


(41) 


(42) 


Finally, equation (40) is used with Equation (42) to give the following 
relation for the EPS mass ratio 



= MRE + MRP + (MRE + MRP + 1 


(1 * HRE) ~ ’) 

I - MRE gf^ 


( 43 ) 


In general, the aV term in Equation (43) comprises all on-orbit func- 
tions such as stationkeeping, repositioning and disposal which are done by 
auxiliary propulsion at the I^p shown in Equation (43). 
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2.3 OrXULTATlON RELATIONS 

Figure 4 shows the circular satellite orbit of radius r around a 
spherical Earth of radius a. The occultation is the fraction of the time f 
that the satellite is in Earth's shadow. 

Figure 4 depicts the Earth and the satellite orbit as seen on the 
plane containing the orbital axis (i.e., the angular momentum vector w) and 
the Earth-Sun line. The orbit viewed toward the Sun is shown in Figure 5. 

The semiminor axis of the perceived ellipse is 

b = r coi 3 (44) 

The equations of the ellipse and the Earth circle in this x-y plane 

are: 


for the ellipse, and 


2 2 


or X + 


2 

cos 3 


= r 


2^2 2 
X + y = a 


(45) 


(46) 


for the circle. 

These intersect at a value of x obtained by simultaneous solution of 
Equations (45) and (46), 


or 

Now, 


cos^B 


n — 2 — T“ 

'a^ - r cos 3 


■ sin 3 
X = r sin uj t = r sin 


2nt 


(47) 
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Figure 4. Plane Containing Orbital Axis and Earth-Sun Line 



Figure 5. View Toward Sun 
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where t is orbit time and T is orbit period. The occultation time is twice 
the time for passage from x = o to x = r sin wt. Hence, the occultation 
fraction f is; 


f 


2t 

T 


1 . -1 
- sin 

TI 



(48) 


The problem then becomes that of findint^ ti, given the orbital para- 
meters and day of the year. It is convenient to adopt the Earth-centered 
coordinate system shown in Figure 6 where the x-y-z cartesian axes are 
referred to the fixed stars, x-y is the ecliptic plane, and z is the 
ecliptic pole. The Earth's axis E is fixed in the x-z plane and the Sun 
line is in the x-y (ecliptic) plane, with angle o equal, in degrees, to 
0.986 times the number of days after the summer solstice. 

The satellite orbit axis w is located by two coordinates; its inclina- 
tion i with the Earth axis and the longitudinal separation of the orbit 
axis and the ecliptic pole. Angle B is the longitudinal separation of E 
and ( 1 ), referred to the ecliptic axis, t = 23.5 degrees, the inclination 
of E. 

The angle P is determined by applying the law of cosines three times 
to the two spherical triangles drawn in Figure 6. Taking P as the angle 
between w and Z, its arc on the sphere is; 

cos P = cos i cos t + sin i sin c cos iji (49) 


Also, since 


cos i = cos t cos P + sin i sin P cos B, 
P as determined in Equation (49) is used to solve for B; 

D cos i - cos e cos P 

COS D - 1 ' ' K 

Sin e sin P 


Finally, from the larger triangle, 

cos n = sin P cos (• + B) 


(50) 

(51) 
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Using t)(deg) = 0.986 times the number of days after summer solstice, 
the orbit inclination i, the inclination of the Earth axis c 23.5 degrees, 
and using iji = -2.1 x 10^^ cos deg/day for the nodal regression rate 

(where r is orbit radius in kilometers), Equations (49), (50) and (51) are 
solved for 3 . The result is used in Equation (48) to calculate the occul- 
tation fraction f. The initial value of ij; depends on the time of day of 
the launch. This is accounted for by directly assigning an initial value 
to this angle. 

The spherical triangle relations used above require that the angles 
be in the range from 0 to 180 degrees. This causes no difficulty in 
Equation (49), with the angles i and c in the first quadrant. Also the 
angle causes no problem. When the computer does not distinguish between 
cos ij; and cos (-ij/), the angle is automatically referred to the proper 
triangle location and geometry for all quadrants, so that P is correctly 
calculated. 

The calculation of B in Equation (50) does require attention. When B 
increases through 180 degrees the computer may interpret the decreasing 
magnitude of the cos as an angle less than 180 degrees. The correction 
procedure, then, is to let B be the value calculated from Equation (50) 
when sin B > 0, and when sin B < 0 the value calculated is replaced by 
(360 deg - B). The program then calculates occultation correctly as the 
climb progress from any starting conditions. 
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2.4 SPACE ENVIRONMENT EFFECTS 


The two space environment effects of importance for orbit transfer 
from Shuttle altitude to geosynchronous orbit are atmospheric drag and the 
Van Allen radiation belt. 

2.4.1 Atmospheric Density 

Figure 7 shows the large variation in atmospheric density that occurs 
at high altitude, even when sunspot activity variations are averaged out 
(Reference 5). For the present study the nominal variation of density 
shown in Figure 7 is used. 

It is seen in this figure that drag effects for transfer to GEO are 
significant only during the early stage of the flight when the altitude is 
below about 400 km. Drag is very important, however, in determining the 
minimum altitude for initiation of the transfer. This subject is dis- 
cussed in Section 2.6. 

2.4.2 Radiation 

The radiation environment is shown in Figure 8 which shows fluence 
rate versus altitude for 0- and 30-degree inclination orbits. Data are 
from Reference 6. The fluence rate is the number of damage events per 
square centimeter per year expressed in terms of the number of 1-mev 
electrons required to produce the same damage. This is a function of cell 
type and cover protection. The data in Figure 8 are for silicon solar 
cells with 6-mil equivalent silica covers, and the fluence rates shown are 
for each side of the cell. For a free standing array, such as would be 
used for present missions, the numbers from Figure 8 are doubled to obtain 
the total damage rate. It is seen that the fluence rate for the equatorial 
orbi^ is more than twice that of the 30-degree orbit. This affects the 
choice of plane change strategy for transfer to GEO. This subject is 
discussed in Section 2.7. 

The effect of the integrated fluence on silicon solar cell power is 
plotted in Figure 9 (Reference 6). 

2.5 SOLAR CELL COVER TRADEOFF 

Tradeoff calculations were made to find the optimum solar cell cover 
thickness for Shuttle orbit to GEO transfer missions. Although this was 
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Figure 8. Fluence Rate vs Altitude for 0 and 30-deg Orbit Inclinations 
(Each side for 6-mil covers - silicon solar cells) 
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Figure 9. Silicon Solar Cell Degradation Caused by Radiation 





done early in the study program before occultation was included in the 
climb program, the conclusions are believed to be applicable to the 
occulted case. 

Table 4 lists the conditions used for the analysis. The solar power 
specific mass is the same function of cover thickness as that in Table 3 
and the range of thruster specific mass brackets the thruster estimates. 

The lowest values of were used with the highest thruster efficiency, 
the highest values with the lowest efficiency and the intermediate speci- 
fic mass was examined with both the highest and lowest efficiency functions. 
In this manner, the full range of anticipated conditions was examined for 
intermediate and extreme cases. These extremes ranged from shortest to long- 
est times at high radiation fluence. 

Results for the four sets of conditions are plotted in Figures 10 
through 13. Thrust time rather than transfer time is shown, reflecting 
the omission of occultation in this calculation. The is not a sensi- 
tive parameter; optimum values are roughly 1500 and 1700 seconds for the 
high efficiency case and 3700 and 3800 seconds for the low efficiency 
function. Cross-plots of time versus cover thickness at these I^p values 

Table4. Range of Parameter (kg/Watt) and n for 
Cover Thickness Tradeoff 


1 

""th 

a 

s 

“e 

Covers 

(mils) 

n Case 

0.005 

0.0050 

0.0100 

1 

III 

0.005 

0.0059 

0.0109 

3 

III 

0.005 

0.0085 

0.0135 

6 

III 

0.005 

0.0250 

0.0300 

12 

III 

0.010 

0.0059 

0.0159 

3 

I III 

0.010 

0.0085 

0.0185 

6 

I III 

0.010 

0.0250 

0.0350 

12 

I III 

0.020 

0.0059 

0.0259 

3 

I 

0.020 

0.0085 

0.0285 

6 

I 

0.020 

0.0250 

0.0450 

12 

I 
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Effect of Ijp and Cover Thickness on Thrust Time at 
« 0.020 Kg/W. Case I n) 
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are shown in Figure 14. Corresponding values of residual power factor 
are shown In Figure 15. 

The 6-mil cover thickness gives thrust times very close to the minimum 
for all combinations of thruster efficiency and specific mass. In addition, 
this thickness results in a residual power factor high enough that further 
improvement would not be worth the additional time. However, the minima in 
Figure 14 are fairly broad so that increased thickness for the less effi- 
cient and heavier thrusters might offer enough improvement in residual 
power to be important for some payloads. In any event, cover thicknesses 
less than 6 mils do not appear to be needed for this orbit transfer 
mission. 

2.6 LOW ALTITUDE FLIGHT 

Conservative calculations of solar array drag effects based on total 
array area, i.e., with the array perpendicular to the flight direction, im- 
pose an unrealistic limitation on minimum useful altitude for solar power. 
The limit occurs well above the Shuttle parking orbit for maximum payload 
(about 250 km). This apparent limitation is readily removed by a feathered 
array strategy in wnich the array is oriented parallel to the flight direc- 
tion and the spacecraft is rolled to the angle for maximin illumination. 

This reduces the available power but the drag area is reduced considerably 
more than the power with the result that the power per unit drag area is 
increased by about a decade. Sustained flight is then possible at alti- 
tudes considerably lower than the Shuttle orbit, and electrically pro- 
pelled spacecraft deployment directly from Shuttle becomes practical. 

2.6.1 Power Factor in Feathered Flight 

The geometric relations defining power factor and roll angle for 
maximum power in feathered flight are developed here. 

Figure 16 shows the geometric relations as seen from the spacecraft, 
i.e., the viewing sphere is centered at the spacecraft. The velocity 
vector pierces the sphere at x which lies on the horizontal plane, the 
zenith is perpendicular to the horizontal plane at high noon and the 
orbit normal pierces the sphere at the right extremity. The angle a is 
zero at high noon and increases to 2v in one revolution; it reaches the 
terminator crossing at n/2 in 1/4 revolution. is shown both as a great 
circle length on the sphere and as a longitudinal angle around the orbit 
normal . 
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Figure 14. Effect of Solar Cell Cover Thickness on Thrust 
Time for Different Combinations of Thruster 
a and n 
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Figure 16. Viewing Sphere Centered at Spacecraft 
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The sun direction makes the angle 90-6 with the orbit plane. The 

sine of this angle is the power factor if the array is not rolled out of 

the orbit plane. 

When the spacecraft is rolled for maximum illumination in feathered 
flight, the normal to the array makes the angle (> relative to the orbit 
plane. This angle is a minimum at high noon and a = 90 degrees (at the 
terminator crossing). When the spacecraft is rolled through the angle ((>, 
the Sun line makes the angle B'.'lative to the plane of the array. 

The spherical right triangle XPS in Figure 16 has sides 90 degrees - 

a, 90 degrees - 6, and hypotenuse e' - Using spherical triangle relations 

it can be shown that 


cos b’ = sin a sin B 


and 


COS a 


(52) 


(53) 


The optimum roll angle is defined by Equation (53) and the power factor 
is sin 6' . 

At high neon, when a = o, cos B' = 0 so that sin b' = 1.0; i.e., the 
array receiver full illumination. At the terminator crossing, a = 

90 degrees and cos B' * sin s. At this location, sin (cos‘^ sin b) ■ 
sin (90 degrees - b) * cos b which is the correct value for the power 
factor when the array is parallel to the orbit plane. 

Referring to Figure 6, the minimum value of b occurs at the summer 
solstice and when ip = 0 degrees. At this point, B = 90 degrees - 
(i + e) “ 38 degrees for an orbit inclination i - 28.5 degrees, b then 
varies within the range from 38 to 90 degrees as the time of the year 
varies and as the orbit regresses. This is shown in Figure 17 where the 
sun angle relative to the orbit plane (90 - fj) is plotted versus the 
number of days after the summer solstice for an assumed altitude of 250 km. 





Figure 18 shows the variation of the power factor (for the optimum 
roll angle 41) as a function of the spacecraft position in orbit, a, and 
for different values of the angle between the orbit normal and the Sun 
direction, for the 250 km example. The value of this power factor, 
averaged over a, is plotted in Figure 19 as a function of The improve- 
ment in power factor at the optimum roll angle is impressive compared to 
cos B which is the power factor for the array parallel to the orbit plane. 
For example at t' = 90 degrees, when the sun is in the orbit plane and 
cos = 0, the average power factor for the optimum roll angle is 0.64. 

Since b varies with time as shown in Figure 17, the unocculted power 
available, averaged over the orbit period, is also a function of time. 

This variation can be calculated by finding b di> a function of time from 
Figure 17 and looking up the corresponding average power factor in Fig- 
ure 19. The result is shown as the upper curve in Figure 20. The occul- 
tation factor, calculated as explained in Section 2.3, is also shown in 
Figure 20 along with the occulted power which is the product of the 
feathered array power factor with the occultation power factor. 

As mentioned earlier, the above illustrations apply for the assumed 
altitude of 250 km. When the equations for occultation and feathered 
flight are incorporated in a programmed calculation of the climb equation, 
the relations shown illustratively are calculated as continuous functions 
of time during the climb. 

It is interesting to continue the 250 km example one step farther to 
estimate the average power available as a function of percentage of the 
time during the year that it is available. Both the occulted and non- 
occulted results are shown in Figure 21. For this example, which is the 
Shuttle parking orbit for maximum payload, 401- of the array power for full 
illumination is available at any time during the year. For 20 % of the 
time, the average power over an orbit period is 50% of maximum power. 

This result is of particular importance for other than orbit transfer 
missions where the spacecraft may remain in a low altitude orbit for a 
substantial time. 
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Figure 18. Feathered Array Power Factor vs Orbit Position 
Angle o for Different Values of 3 (no 
occultation) 
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Figure 19. Average Power Factor vs (i for Feathered Array 
(No occulation) 
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Figure 20. Average Power Factor vs Time After Summer Solstice 
for Feathered Array (28.5® Orbit Inclination, 

H » 250 km) 
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Figure 21. Percent Time Available vs Average Power Factor Per 
Day for Feathered Array (28.5° Orbit Inclination, 

H = 250 km) 
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2.6.2 Orbit Transfer Starting Power 

For the orbit transfer mission, the spacecraft climbs quickly out of 
the drag region so that power loss associated with orbit regression is not 
an important consideration. 

Referring to Figure 6, the optimum time of day to start the transfer 
is that for which B = -0 so that the orbit vector is coplanar with the 

ecliptic pole and the sun direction, as illustrated in Figure 22. 

Equations (64) and (55) result from the law of cosines applied to the 
small triangle in Figure 22. 

cos F = cos i cos e + sin i sin r cos <|> (54) 

cos i = cos e cos P + sin e sin P cos e (55) 

Since P = 90 - B, 


cos P * sin B 

sin P = cos B 

Equations (54) and (55) can be rewritten as: 

sin B = cos i cos e + sin i sin e cos ^ (56) 

ft _ cos i - cos’ € sin0 

Then with i and e given, <|f is assigned to determine B from Equation (56) 
and this is used in Equation (57) to calculate 8. In this indirect way, 

B is found as a function of 0 and Figure 19 is used to find the corres* 
ponding power factor averaged over the orbit. The occultation effect is 
accounted for in the ascent calculation. 

Figure 23 shows the result of this calculation plotted as maximum 
average power factor at the start of the orbit transfer versus time in 
days after the sunnier solstice. The unocculted power available from the 


48 


mm 






OF ^ ' 





ORIGINAL PAGE IS 
OF POOR QUALITY 


feathered array varies from 70% to 90% of the full array power. For 2/3 of 
the year the available power is more than 80% of the full array power. 

2.6.3 Sustained Low Altitude Flight 

Both Figures 21 and 23 show that the feathered array, rolled for maxi- 
mum illumination during an orbit, produces an average power of about 70% of 
the full array power. The starting power for the transfer mission is 
about 10% higher than the average value for sustained low altitude flight. 

Figure 24 shows that the occultation fraction is not very sensitive 
to altitude in the neighborhood of 250 km where drag effects are large 
enough to require feathering the array, i.e., from 100 to 400 km. There- 
fore, a reasonable assumption is that an average power over the orbit 
period equal to 40% of the full array power is available over this alti- 
tude range and at any time of the year, for either low altitude sustain<>d 
flight or for the low altitude end of an orbit transfer mission. The cor- 
responding drag area in feathered flight is assumed to be 5% of the total 
array area. 

The power per unit drag area is therefore increased by a factor of 8, 
from 150 to 1200 kg/W, in the feathered configuration as compared to the 
full blanket drag area and normal incidence power. This not only reduces 
the drag effect during the early stage of an orbit transfer to GEO, but it 
greatly reduces the minimum altitude for sustained flight. 

This minimum is found by setting the numerator in Equation (16) equal 
to zero. After some rearrangement, the following expressior' is obtained 
for the required electric propulsion mass ratio: 


PL 




2 r \ 


2 - 



1 


( 58 ) 
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The drag contribution of the payload area appears in the numerator 
ratio which shows that the electric propulsion mass ratio increases 
linearly with the power supply specific mass and inversely with the r»»y- 
load specific mass. The drag effect of the power supply area appears in 
the denominator where it determines the condition for which the denomina- 
tor becomes zero and the electric propulsion mass for equilibrium flight 
becomes infinite. 

Equation (58) was used to plot the minimum altitude relations shown 
in Figure 25 for nuclear and solar power. The nuclear power source was 
estimated to have a specific mass = 0.033 while for solar power 
ranged from 0.012 to 0.036 (kg/W). 

The most important conclusion to be drawn from Figure 25 is that 
sustained flight with electric propulsion is practical down to 150 km, 
i.e., 100 km below the Shuttle parking orbit for maximum payload delivery. 

In this high drag region, the feathered solar array performs as well ai 
the nuclear source. Thus, solar electric propulsion can be usod directly 
to augment Space Transportation System capability without relying on inter- 
mediate chemical propulsion stages. Solar electric orbit transfer vehicles, 
servicing vehicles, or on-board propulsion systems may be deployed directly 
from the Shuttle Orbiter. 

2.7 PLANE CHANGE 

Edelbaum (Reference 7) has shown that the characteristic velocity aV 
required for low thrust orbit transfer with continuous plane change i, at 
constant I^p, can be expressed by the equation, 

aV * - 2vjVj. Cos J i + Vp (59) 

where Vj and Vp are initial and final orbit velocities defined by Equa- 
tion (1). By combining Equations (1) and (59), the characteristic velocity 
is expressed dimensionlessly in terms of the initial and final orbit radii 
and plane change angle, 
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This equation can be used to compare four different ways to change the 
orbit plane angle during transfer from a low earth orbit of 250 km to GEO, 
as shown in Figure 26. The upper two curves compare the very inefficient 
approach of making the plane change at the initial altitude, with the more 
efficient one in which the plane change is made at GEO. The smallest aV 
requirement occurs when the plane rotation is done continuously during the 
transfer, according to an optimum thrust vector orientation schedule given 
in Reference 7. Although this minimizes AV, it has the disadvantage of 
producing a greater total radiation fluence for the case of a 28.5 degree 
launch followed by a plane rotation to 0 degrees. This effect is seen by 
comparing the fluence rate curves in Figure 8 for 0 and 30 degree inclina- 
tion. The ratio of these rates increases from 1.0 at 1500 km to 2.3 at 
maximum fluence rate (5000 km) and to 2.9 at 10,000 km. 

This suggests the alternate plane change approach in which the orbit 
inclination remains fixed up to 10,000 km after which it is rotated to 
the final value according to the optimum schedule. Figure 26 shows that 
the aV penalty relative to the optimum plane change is small. It is about 
1/3 of that due to making the total plane change at GEO. In terms of Av, 
the transfer from 250 km to GEO, with i changing from 28.5 degrees at 
launch to 0 degrees at GEO, requires 5990 m/sec for the o|. l:imum schedule 
while an additional 349 m/sec (5.8%) more aV is needed when all of the 
change is made above 10,000 km. 

R. DiEsposti, at NASA-LeRC, has calculated the power loss for both the 

optimum steering program and the alternate approach with plane change above 

10,000 km. For a 3-mil shield thickness with infinite back-shielding, he 

found the residual power advantage of the alternate approach to be less than 

★ 

2% greater than for optimum steering. This suggests that optimum steering 
may be the preferred approach, although detailed trade calculations would 
be needed to verify this for a specific spacecraft design. 


* 

Private communication. 
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2.8 ORBIT TRANSFER TIMES 


Orbit transfer times presented here include all propulsion 
requirements for transfer from LEO to GEO followed by 5 and 10 year life- 
times and removal to a disposal altitude. Effects of occultation, drag and 
solar cell degradation are included. On orbit aV requirements are discussed 
in Section 3. 

2.8.1 Nominal Mission Performance 

The effect of thruster efficiency and I^p on orbit transfer time is 
analyzed for different constant values of the electric propulsion system 
mass ratio. The calculation includes the effects of atmospheric drag, solar 
cell degration and occultation for mission parameter sets shown in Table 5. 


Table 5. Mission Parameter Sets 


Mission Set No. 

/"eps\ 

Wl/t 

Lifetime on Orbit 
(years) 

aV on Orbit 
\sec/ 

1 

1.0 

10 

1731 

2 

2.0 

10 

1731 

3 

3.0 

10 

1731 

4 

4.0 

10 

1731 

5 

2,0 

5 

949 


Initial altitude, -250 km 

Final operation altitude, = 35786 km 

Initial inclination, i^ = 28.5° 

Final operation inclination, i^ = 0° (equatorial orbit) 

Disposal altitude, = 40786 km 

( ^EPS I represents the ratio of the total propulsion mass (primary 

^PL /y thruster system dry, primary propellant and auxiliary propul- 
sion required for stationkeeping and disposal) to the 
payload mass. 
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The spacecraft parameters are; 

= 0.024 kg/W (thruster, solar arsay and processor) 

2 

= 150 W/m (solar array) 

^PL 2 

A — = 100 kg/m^ (payload) 

MpL 


n = 20 to 100% 


'SP 


M 

Vl 


= minimum to 5000 seconds 

Mpp 

and Tj — , electric propulsion dry mass ratios as required for the 

ripL 

total electric propulsion system mass ratios in 
Table 5. 


Figures 27 through 30 show the effect of thruster efficiency and spe- 
cific impulse on orbit transfer time for total electric propulsion masses 
of 1, 2, 3 and 4 times the payload mass for 10-year missions. This mass 
accomplishes the total maneuver requirement of orbit transfer, inclination 
change, attitude control, stationkeeping and disposal. Each transfer time 
curve on these figures exhibits an optimum. To the left of the optimum, 
propellant mass requirements at low specific impulse are most influential. 
To the right of the optimum, power requirements at high specific impulse 
are more significant. 


Figures 27 through 30 show that an optimum I^p oi b that produces 
the shortest trip time for a given thruster efficiency; as seen in the 
figures, it is independent of t\. Also plotted in these figures is the 
Case II efficiency relation, which shows the characteristic drooping 
toward lower I^p's. The net effect of this droop, as shown in Table 6, 
is that minimum trip time occurs at an higher than the "optimum" value. 

A plot of Ijp and corresponding minimum transfer times is shown in 
Figure 31 as a function of payload fraction, i.e., the ratio of payload 
mass to total initial spacecraft mass. This fraction is related to the 
total electric propulsion mass ratio as shown below: 


Payload Fraction = 



+ 1 


58 









61 








OR!G[:;;.L t'. / - 

OF POOR Q 


Y 



Figure 31. Effect of Payload Fraction on Minimum Transfer Time and 
Corresponding I^p (Case II n, a = 0.024 kg/W 10 Years) 
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Table 6. Comparison of Optimum I^p and I^p For Minimum 
Transfer Time - Case II Efficiency 


/”eps\ 

Wl 

Optimum I^p 
(sec) 

Ijp For Minimum 

Transfer Time 
(sec) 

1 

2600 

3100 

2 

1750 

2300 

3 

1550 

2000 

4 

1400 

1800 


It is shown in the figure as a percentage rather than a ratio. 

From a system point of view, it is more economical to operate at an 
Ijp slightly lower than the minimum transfer time value in order to 
reduce the more expensive electric propulsion mass while increasing the 
less expensive propellant mass. 

The relative magnitudes of these mass tradeoffs are seen in Fig- 
ures 32 through 35 for the 10-year mission cases shown in Figures 27 
through 30. The preferred I^p conditions are shown as minimum transfer 
times in the figures. 

The transfer time results for Mission Set 5 are plotted in Figure 36 
and the mass ratio breakdown is shown in Figure 37. The greatly reduced 
on-orbit propellant for the shorter on-orbit life is apprent. It is 
also seen in Figure 38, which also shows the effect of payload fraction 
on the percent of total propellant that is used for on-orbit functions. 

The reduction of on-orbit propulsion from 12.6 to 7.3% of the total 
propellant mass results in a small reduction in transfer time. As shown 
in Figure 39, the transfer time for the 5-year mission is about 6% less 
than for the 10-year mission. 

The two thruster parameters that have the greatest impact on transfer 
time are thruster efficiency and the electric propulsion specific mass 
parameter a. Figure 40 is a cross-plot of Figures 27 through 30 showing 
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Figure 33. Propulsion Mass Breakdown for 25t Payload Fraction 
(Case II p,, a = 0.024 kg/W, 10 Years) 
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Figure 35. Propulsion Mass Breakdown for 50% Payload Fraction 
(Case II n, a = 0.024 kg/W, 10 Years) 
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Figure 37. Propulsion Mass Breakdown for 33% Payload Fraction 
(Case II n, o( = 0.024, 5 Years) 
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Figure 38. Percent of Total Propellant Used On-Orbit 
(Case II n. o = 0.024, I^p for Minimum Transfer Time) 
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transfer time as a function of thruster efficiency at the optimum I^p for 
the four different values of the total electric propulsion mass ratio. The 
corresponding optimum I^p values are also shown. For the case of mass 
ratio equal to 2.0 at 50% efficiency, for example, the slope corresponds 
to a 3.7-day reduction in transfer time for a 1% improvement in thruster 
efficiency. 

Figure 41 shows that transfer time varies essentially linearly with 
a. For a total electric propulsion mass ratio of 2.0 and a thruster 
efficiency of 45% at 1750 seconds I^p, the slope at a = 0.024 kg/W 
indicates a ti ansfer time reduction of 20 days for a 10% reduction in a. 

t 

Figure 42 shows the power remaining after orbit transfer as a function 
of total EPS mass ratio, assuming the climb is conducted at optimum I^p and 
Case II efficiency. In the mass ratio range of 2 to 4, the residual 

power ranges from 55 to 59% of the beginning of life value. At a mass ? 

ratio of 1.0, it is 50% of tlie initial value. The on-orbit power degrade- ■ 

tion is small. Table 7 compares the integrated fluence over a 10-year 

mission with that accumulated during orbit transfer and with tne total 

fluence per 11 year solar cycle. The total on-orbit increment is about 

2% of the orbit transfer value. Solar cell covers equivalent to 6 mils 

of fused silica were assumed for this calculation. 



Table 7. Comparison of Orbit Transfer and On-Orbit Fluence Levels 
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Figure 41. Effect of Specific Mass on Transfer Time for 
Case II n = 45% at 1750 sec 
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Figure 42. Residual On-Orbit Power for Ootimum I^p 
Transfer and Case II Efficiency 

(o = 0.024 kg/W, p = 150 W/m^, Mp^/Apj^ = 100 kg/tn^, 10 years) 
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2.8.2 Fast Orbit Transfer 


The great attraction of electric propulsion is its ability to 
transfer large payload fractions and its greatest drawback is long transfer 
time. It is interesting to examine the: payload penalty associated with 
fast orbit transfer for advanced technology thruster having the low 
values of a and the improved efficiency that might be anticipated. 

The resistojet offers high efficiencies, in the range of 80 to 90% 
but is limited in I^p to values less than 1000 seconds. The pulsed induc- 
tive thruster operates in the neighborhood of 50% efficiency in the range 
of 1500 to 2500 seconds. Both of these thruster types have lightweight 
power conditioning and appear to be possible candidates for fast transfer 
from Shuttle orbit to GEO. 

To assess this possibility quantitatively, the parameter a was reduced 
from the value 0.024 kg/W, used previously, to 0.012 kg/W, the Case III 
efficiency function was used for the pulsed inductive thruster and two 
resistojet cases were examined: 90% at '900 seconds I^p and 100% at 

1000 seconds I^p. The payload fraction was varied frcxn 10% to 50%. 

Figures 43 through 48 show the results of the transfer time calcula- 
tions, assuming a 10-year duration on-orbit; the corresponding propellant 
breakdown plots are shown in Figures 49 through 54. These figures have the 
same characteristics seen in the previous section. Note, however, the steep 
slope of the transfer time curves at low specific impulse. Mission sensi- 
tivity to thruster operating characteristics in this region is extremely 
high. This will require strict adherence to performance specifications for 
any thruster operating in this range. 

Figure 55 summarizes the 10-year on-orbit propellant fraction used for 
stationkeeping and disposal. The values range from 6% of the total pro- 
pellant at 10% payload fraction, to 16% at 50% payload fraction. This 
increase is due to the decreased orbit transfer propellant requirement 
that results from the increase in I^p for minimiR transfer time as the 
payload fraction increases. This trend is seen in Figure 56 which shows 
the effect of payload fraction on the minimum transfer time and also on 
the corresponding I^p for the assumed Case III efficiency. 

At 10% payload fraction the l^p for minimum transfer time is 1300 sec- 
onds and the time is only 31 days. This condition, however, is unrealistic 
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Figure 46. Transfer Time Trade-Off - 301C Payloed Fraction 
(a « 0.012 kg/U. 10 Years) 
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Figure 49. Propulsion Mass Breakdown for 10% Payload Fraction 
(Case III n, o = 0.012 kg/W, 10 Years) 
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Figure 51. Propulsion Mass Breakdown for 20% Payload Fraction 
(Case III n, a » 0.012 kg/W, 10 Years) 
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Figure 53 Propul <;1on Mass Breakdown for 40% Payload Fraction 
(Cas? IIT r,, a = 0.012 kg/W, 10 Years) 
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Figure 54. Propulsion Mass Breakdown for 50% Payload Fraction 
(Case III n» o = 0.012 kg/W, 10 Years) 
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Figure 55. Percent of Total Propellant Used On-Orbit 
(Case III n. o - 0.012 kg/W. 10 Years) 
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because payload fractions as low as 10% can be managed by chemical propul- 
sion with transfer times measured in hours rather than days. To be attrac- 
tive, electric propulsion must transfer payload fractions that are 
significantly greater than the maximum capability of chemical propulsion. 

A reasonable lower limit for payload fractions of interest for 
electric propulsion is 20%. Figure 56 shows a 40-day transfer time require- 
ment for the plasma thruster at this payload fraction. Transfer time 
Increases rapidly with payload fraction; at 30% it is 52 days, it is 75 days 
at 40% and at 50% payload fraction it has increased to 123 days. However, 
even these longer times are short compared to those previously indicated 
for conventional electric propulsion. The transfer time reduction of 
more than 50% is due both to the high thruster efficiencies of the 
Case III efficiency relation and the low value of specific mass a for 
the pulsed inductive plasma thruster. 

The resistojet is capable of high efficiencies at Ijp values that are 
large compared to chemical rockets and low relative to electric thrusters 
that use either magnetic or electric fields for propellant acceleration. 
Conversion of thermal to kinetic energy is of the order of 90% until the 
Ijp limit fixed by material properties and heat loss is approached. For 
hydrogen propellant this might be about 900 seconds. Thus the resistojet 
range for the performance maps in Figures 43 through 48 is co the left of 
the optimum I^p. 

Figure 57 summarizes these zones of the performance maps for transfer 
times of 20, 30, 50 and 100 days. Again, 20% payload fraction is taken to 
be the minimum value of Interest for electric propulsion. This limit and 
the maximum hp of 900 seconds are shown In Figure 57 for 90% thruster 
efficiency. For 20% propellant fraction the trip time is about 30 days. 

If the Ijp were reduced to 800 seconds, the time would be 50 days. 

It Is also seen In Figure 57 that trip time Increases more rapidly 
with payload fraction for the resistojet than It does for the pulsed 
1ndix:t1ve thruster. For example, at 30% payload fraction it has increased 
to 100 days as compared to 52 days for the pulsed inductive thruster. 

Figure 58 shows this effect as a function of payload fraction for 90% 
efficiency at 900 seconds and also for the limiting case of 100% efficiency 
at 1000 seconds. 
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Thus with forseeable increases in thruster efficiency and reduction 
in specific mass (Case III efficiency relation and a ® 0.012 kg/W), Shuttle 
orbit to GEO transfer time can be reduced to 40 days for 20% payload 
fraction. At 50% payload fraction, the corresponding transfer tin»o is 
only 123 days. Resistojet performance is only attractive for LEO to GEO 
transfer missions at low payload fractions. Us transfer times quickly 
become prohibitive with increasing payload fraction. 
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3. INTEGRATED PROPULSION ANALYSIS 


The purpose of this investigation is to identify potential advantages 
of using electric propulsion not only for orbit raising to geosynchronous 
orbit but also for additional (auxiliary) functions that will be required 
during the ascent and orbital phases of the mission. These functions can 
be accomplished by using the EPS as an integrated part of the space vehicle 
rather than as a separate orbit transfer vehicle only. 

3,1 SCOPE OF STUDY 

For the integrated propulsion system, functions to be considered in 
addition to orbit transfer include: 

1) Stationkeeping 

2) Orbital repositioning 

3) Spacecraft disposal at the end of the mission to an orbital 
altitude that precludes future interference with other satellites 

4) Attitude control, including correction of perturbing torques due 
to gravity gradient and aerodynamic effects 

5) Propulsion power sharing W'th the payload and with spacecraft 
housekeeping subsystems 

Of these, the first three functions involve translation aV maneuvers, the 
fourth function, rotational maneuvers. 

Results of the analysis are expressed in terms of weight and cost 
savings that result primarily from reduced propellant expenditure on 
orbit with the concommita.it, even greater, reduction of propeilant mass 
required for LEO to GEO transfer. It is these savings that translate into 
a reduction of total launch weight and, consequently, in Shuttle transporta- 
tion cost savings. 

Another important aspect of integrated propulsion is the question of 
thruster size and power requirements for adequately performing the 
auxiliary thrust functions, especially the correction of large distuitance 
torques. Also, the impact of sharing available power effectively between 
thrust operations and operation of other power consumers must be 
co.isidered. 
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In addition to the integrated electric system, the use of a chemical 
auxiliary propulsion system is considered and the possible combination 
of elect»*ic and chemical thruster operations is also examined. 

3.1.1 System Concept and Assumptions 

The principal considerations in this investigation of integrated EPS 
functions are the following: 

1) The orbit transfer from LEO to GEO is performed by means of a 
primary electric propulsion system. 

2) Of primary interest are spacecraft of a size, i.e., total mass 
and dimensions, consistent with a single Shuttle launch rather 
tha:t requiring multiple launches and assembly in orbit. The 
results, however, should be applicable for spacecraft requiring 
several launches. 

3) The primary electric propulsion system is designed to exercise 
thrust vector control, by thruster gimballing or thrust vector 
steering, as well as three-axis attitude control that can be 
done with moderate control torques. Auxiliary control will be 
provided to correct large disturbance torques during primary 
thrust operation and, in general, to perform attitude control 
during primary thrust off-periods. 

4) Questions pertaining to attitude control of specific spacecraft 
configurations are generally not addressed. Also, questions of 
large structure deformation due to thrust forces and torques, and 
dynamic interactions, addressed in other recent studies, e.g. , 

by Boeing Aerospace Company, are considered beyond the scope of 
this study task (Reference 2). 

5) It is assumed that a major portion of the power available 
initially during LEO to GEO transfer will still be available 
for use by the EPS during the orbital mission phase, i.e., 
solar array degradation due to radiation damage during transfer 
will be about 4u%. 

6) Payload and spacecraft housekeeping functioris may use a major 
part of the power available on orbit, thus imposing possible 
constraints on EPS thrust operations by power sharing. 

7) North-south and east-west stationkeeping are to be performed 
without causing a change in the spacecraft nominal orientation, 
while orbital repositioning can be done by the main thrusters 
utilized during orbital transfer, if necessary by reorientating 
the spacecraft from its normal stationkeeping attitude. Final 
disposal can be done by either the main thrusters or the 
auxiliary stationxeeping thrusters. 
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8) Primary as well as auxiliary thrusters must be located and 
orientated so as to avoid or minimize thrust plume impingement 
on spacecraft appendages such as solar panels. However, 
details of thruster placement and articulation are not 
addressed during this task, 

9) Attitude control torque requirements to control major 
disturbances are derived by extrapolation from existing data 
in the literature, e.g.. Reference 2. 

10) The state of the art Isp for auxiliary electric propulsion, 

3000 seconds, is used for all calculations. Optimization to 
other values is not useful because of the low ratio of 
secondary to primary propulsion total impulse. 

3.1.2 S pacecraft Description 

This study makes use of results obtained from previous and concurrent 
studies of auxiliary electric propulsion systems on large space structures 
performed by Boeing Aerospace Company (Reference 2) under NASA Lewis 
Research Center Contract. Recent studies by TRW under NASA-LeRC and USAF 
Rocket Propulsion Laboratory contracts of auxiliary propulsion systems 
(References 8 and 9) also provide useful source data. 

Data from the Boeing APS study used in this task include: 

t Mass properties of various spacecraft structures 

• Disturbance torque characteristics 

t Attitude control and maneuver requirements 

• Thruster sizing and placement 

• Propulsion system design characteristics and propellant 
mass estimates 

A wide range of spacecraft configuration categories, masses and 
dimensions are covered in the Boeing report, with emphasis on large 
structures that require multiple Shuttle launches and assembly in orbit. 
For the present analysis data applicable to a class of Intermediate-sized 
vehicles with nearer-term mission application prospects were selected. 

The data used pertain to a spacecraft category termed "modular 
antenna structure" (see Figure 59) with a mass ranging from 2000 to 
27,000 kg. Spacecraft ch iracteristics in the lower mass range are 
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compatible with a single Shuttle launch. The one-axis gimballed solar 
arrays provide primary propulsion power of several hundred kW. This implies 
array sizes of several thousand square meters with tip-to-tip dimensions 
of hundreds of meters. Data for this spacecraft class are presented in 
the Boeing report with antenna size as scaling parameter (d ranging from 
15 to ?00 m). Corresponding solar array dimensions range to 70 m panel 
length. Results obtained by Boeing are assumed to apply also to other 
intermediate-sized spacecraft of similar configuration where the payload 
is not necessarily a communication system or space radar. 

The Boeing study concludes that the propellant mass required for 
auxiliary propulsion (E-W and N-S stationkeeping, solar pressure compensa- 
tion and attitude control maneuvers) is determined primarily by the aV 
maneuver requirements, with attitude control expenditures being a small 
percentage of the total.* In the present study the aV maneuver require- 
ments are increased by about 50% to include orbital repositioning and 
disposal maneuvers so that the relative proportion of attitude control 
propellant mass is thus further decreased. Consequently, attitude control 
propellant mass was neglected in the total propellant mass estimates. 

Results of earlier TRW auxiliary propulsion integration studies 
(References 10 and 11) of existing satellite configurations designed for 
operation in geosynchronous (Intelsat V, MMS, DSCS III, DSP) and inter- 
mediate orbital altitudes (GPS) are applicable to this study in the 
following areas: 

• Thruscer location and operational sequence 

f Thrust plume impingement avoidance 

• Attitude control propellant mass data 

t EMI and RFI interactions 

These d'ta are useful in auxiliary thruster placement and orientation 
trades, ’n bracketing propellant mass requirements and in defining areas 
of 'urcher study in controlling potential propulsion system/spacecraft 
interactions. 


*Pr1vate communication from W. Smith, Boeing Study Manager. 



3.2 PROPULSION SYSTEM DEFINITION 

The integrated propulsion system comprises the main thrusters used 
for orbit transfer and an assembly of small thrusters to provide the 
much smaller forces needed for stationkeeping and the torques required 
for attitude control. Several configurations were examined before 
selection of the system described here. 

3.2.1 Thruster Arrangement 

Figure 60 shows the locations and orientations of the various 
electric thrusters used or the generic spacecraft design analyzed in 
this study. 

A set of gimballed main thrusters is located at one end of the 
spacecraft center body, with the thrust axis normally pointing along the 
length axis (z-axis), such that propellant depletion does not produce a 
major thrust vector misalignment. 

Four auxiliary electric th»*usters are grouped around the main 
thrusters, pointing normal to the centerline along the x and y axes. 

Two additional clusters of f've thrusters each are located on out- 
rigger booms extending in +x and -x directions with thrust orientations 
as indicated in the figure. The booms must be sufficiently long to pre- 
vent direct thrust plume impingement on nearby portions of the solar panels 
when these panels are oriented parallel or nearly parallel to the space- 
craft x-y plane. The locations shown also provide favorable attitude 
control moments. 

The X, y and z body axis designations shown in the sketch conform 
with the notation used in the Boeing study. The thrust designations 

indicate in which direction the thrust forces (F), 

A y y z 

and around which axes the thruster torques (T), are acting. These desig- 
nations, too, conform with the Boeing study. 

The auxiliary tr.»“usters shown can be operated in sets to produce 
ei\i',er pure translation .nlong any of the body axes or pure moments 
“ground the axes. For example, if pure translation maneuvers without 
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rotation are required the thrust magnitude or pulse length of thrusters 
used in combination must be controlled such that opposing torques cancel 
out. If pure rotation around a principal body axis is desired, care must 
be taken to cancel coupling effects by appropriate combination of several 
thrusters. Rotations around the x, y and z axes are referred to as pitch, 
yaw and roll, respectively. The thruster set is arranged to provide 
partial redundancy. 

Gimballing of the main thrusters can correct thrust axis misalignment 
and can also be used to generate three-axis control moments. For pitch 
or yaw control the thrusters are gimballed jointly either around the x or 
the z axis. Roll control requires differential gimballing of at least 
two main thrusters. The auxiliary thrusters are assumed to be rigidly 
mounted at their respective locations. 

Generally, during the primary thrust phases, most attitude control 
functions are performed by the main thrusters without requiring assistance 
from the electric APS thrusters. However, the control torques generated 
by the primary or auxiliary electric thrusters will not always be large 
enough to overcome maximum external disturbing torques such as may arise 
from gravity gradient or aerodynamic effects at low altitude. These 
conditions will occur for spacecraft with large appendages such as the 
solar arrays with tip-to-tip dimensions typical for the generic spacecraft 
configuration assumed here. Such spacecraft are characterized by large 
differences in maximum and minimum moments of inertia and large center- 
of-pressure/center-of-mass offsets. The moment of inertia differences 
are multiplying factors in the expressions for gravity gradient torques 
while the c.p./c.m. offset governs aerodynamic disturbing torques. 

The Boeing study has shown that auxiliary chemical thrusters may 
have to be added to provide the necessary corrective torque capability 
which in the low and intermediate altitude region can be several orders 
of magnitude greater than that achievable by electric auxiliary thrusters 
of reasonable size. Data to substantiate this requirement for spacecraft 
mass and size characteristics considered here are shown in the next 
section. 
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Another factor to be considered is the periodic power cutoff due to 
eclipses that occur at low as well as high orbital altitudes. Although 
momentum storage devices can maintain normal attitude control functions 
during eclipse intervals, the use of chemical backup thrusters can 
effectively augment this capability. The alternative of using batteries 
with storage capacities of the order of 10 to 30 kW-hr to support 
auxiliary thrust operation does not appear to be practical or cost 
effective (see discussion in the next section). 

Placement of the chemical thrusters at the locations shown in 
Figure 60 serves to control the largest external disturbance torques, 
namely those around the x and z axes. A minimum of six chemical thrusters 
is required. 

3.2.2 Thrust Modes 

In some instances during the ascent and the orbital mission phases, 
the desired maneuver modes require the thrust to be aligned with the 
orbital velocity vector; in others, an out-of-plane orientation at some 
angle or even normal to the velocity vector is required. 

3. 2. 2.1 Ascent Phase Thrust Orientation 

During the ascent phase the main thrusters operate continuously 
except during eclipses. The initial part of the transfer takes place 
without plane change to minimize exposure to Van Allen belt radiation. 
During this part of the ascent phase the thrust vector is oriented along 
the velocity vector. 

During the final two-thirds of the ascent the thrust vector is 
oriented at an increasing out-of-plane angle to accomplish the required 
plane change into equatorial orbit by the time the vehicle reaches 
synchronous altitude, while the principal in-plane thrust component 
continues pointing along the velocity vector. Periodic reversals of the 
out-of-plane orientation are performed twice per orbital revolution, 
according to a steering program that places the midpoints of each succes- 
sive out-of-plane deflection at tne ascending and descending nodes of 
the orbit (see Reference 2, Section 4). In the preferred ascent profile 
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the rotational rates and accelerations around the local vertical that 
are required to achieve this steering program, however, are only about 
25% of the maximum values quoted in Reference 2, since the out-of-plane 
maneuvers do not start before the orbital period has increased to about 
6 hours. The required lateral steering torques can be provided by main 
thrust gimbal deflection assisted if necessary by the appropriate 
auxiliary thrusters. 

In the low altitude region of ascent, periodic roll maneuvers are 
required to maintain optimum sun exposure of the feathered solar array 
panels (Section 2.6.1). This may present problems of excessive roll 
torque requirements since the roll moment of inertia is the largest 
one for the reference configuration being considered, and the maximum 
roll rates required to track the sun and achieve roll excursions of 
±75 degrees twice per orbit are as high as (±) 0.03 deg/sec. However, 
this occurs only under worst-case conditions, i.e. , when the sun is within 
about 15 degrees of the orbit plane. The problem may be avoided by 
restricting the launch date and time-of-day so as to get a more favorable 
sun alignment, and/or by limiting the roll excursion at the expense of a 
small power loss. 

A performance tradeoff between roll maneuver requirements for optimum 
sun exposure versus degraded power generation, and also versus required 
propellant expenditures for chemical thruster operation to support these 
maneuvers would be needed to define the best operation strategy. 

3. 2. 2. 2 On-Orbit Maneuvers 

On-orbit operation normally requires the spacecraft z axis (length 
axis) to be Earth pointing, e.g., to perform typical Earth observation or 
communicatioi.a relay tasks. This constraint prevents the use of the main 
thrusters for nonradlally oriented aV maneuvers, such as E-W or N-S 
stationkeeping. 

Figure 61 shows the spacecraft orientation during the orbital phase 
with the z axis pointing along the radius vector to earth, the x-axis is 
In the oi^bit plane. This orientation is consistent with best sun orienta- 
tion of the solar panels which must perform one full revolution around 
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their center line (the y-axis) per day. It is used by most three-axis 
controlled synchronous equatorial satellites today. 

It is seen from Figure 60 that N-S stationkeeping, which requires a 
aV maneuver normal to the orbit plane, is performed by a combination of 
Fy thrusters, generally in sets of three, to nullify any pitch or roll 
torques that would be induced by individual thrusters. E-W stationkeeping 
requires thrust along the velocity vector or in the opposite direction 
and is performed by firing F^ thrusters. 

The solar panels are in a worst case orientation for plume impinge- 
ment from the boom-mounted F^ thrusters once every 12 hours, and s1mila»‘ly, 
but 6 hours later, for plune impingement from the F^ \,hrusters mounted 
near the main thruster assen>bly. These thrusters thus may require mount- 
ing at a fixed outward cant angle, depending on their position relative to 
the solar panel edges. 

For satellites with large Earth pointing antennas, centered on the 
z-axis, a similar outward cant angle may be required for the forward firing 
F^ thrusters. 

The required cant angle varies inversely with boom length. A trade 
between extra boom structure mass and extra propellant mass expenditure 
due to the cant angle is involved in determining the best arrangement for 
plume impingement avoidance. 

Orbital repositioning maneuvers to a different longitude may be 
required several times during the spacecraft life. For these maneuvers 
the major aV Impulses, one In posigrade and one In retrograde direction, 
are required at the beginning and the end of the transfer In a sequence 
that depends on whether eastward or westward respositloning Is Intended. 
Since normal payload operations, e.g., communications relay functions, 
can be suspended during the repositioning phase which may last for 
several weeks, the vehicle can be reoriented again to permit utilization 
of the main thrusters. 

3.2.3 Thruster Function Allocation 

From the foregoing discussion It Is apparent that the use of 
auxiliary electric thrusters Is a practical approach to performing 
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translational and rotational maneuvers during the orbital mission phase 
and to providing attitude control support during the ascent phase. How- 
ever, chemical backup thrusters may be needed; (1) to meet peak control 
torque requirements that exceed the capability of auxiliary electric 
thrusters of reasonable size, and (2) to substitute for electrical 
thrusters under eclipse conditions (unless momentum storage devices can 
fill this need). 

Table 8 lists the thruster types required and their utilization modes 
versus mission operation phases, to summarize the results of this discus- 
sion. Table 9 summarizes the thruster locations and orientations assumed, 
their functions in translational and rotational maneuvers and the minimun 
number of thrusters employed. 

3.3 SYSTEM REQUIREMENTS 

This section presents performance data, design characteristics and 
other parameters needed in assessing the applicability of integrated 
auxiliary electric propulsion as compared with an alternative chemical 
APS. These data are required to substantiate assessments previously 
quoted from the Boeing APS study (see Section 3.1) and also, to support 
propellant mass evaluations and related analyses in Section 3.4. Appli- 
cable data derived from the reference study are presented as appropriate. 

3.3.1 Spacecraft Mass Characteristics 

Table 10 lists the net spacecraft mass and moments of Inertia of the 
selected reference spacecraft class ("modular antenna system") and two 
other spacecraft classes of larger and smaller size investigated in the 
Boeing study. These spacecraft categories are defined by Boeing as "large 
erectable antenna structures" and "concentrated mass with appendages,” 
respectively. In each case the mass and moment-of- inert la values are 
listed for the parametric range investigated. The mass characteristics 
of the selected reference system fall into an Intermediate range between 
those of the other two. 

Figure 62 shows spacecraft masses and the largest and smallest 
moments of inertia, I^^ and 1^^, versus antenna diameter, d, in the range 
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Table 8. Thruster Function Allocations 
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Table 9. Auxiliary Thruster Functions and Numbers 
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Table 10. Mass Properties Data for Three Spacecraft/Structure Classes 

(Reference 2) 














vr 


JO ry^p<4^0AV 


111 


Figure 62. Mass and Moments of Inertia of Reference Spacecraft/Structure (Modular 
Antenna Structure) Extrapolated from Boeing Oata (Reference 2) 







ot lb to 200 m. The scaling relations used to express other spacecraft 
characteristics in terms of antenna size are listed in Table 11. The an- 
tennas were taken to be solid, as opposed to mesh, construction for the fol- 
lowing analysis. 

3.3.2 Disturbance Forces and Torques 

Figures 63 through 65 show disturbance forces in the x, y and z direc- 
•.ions acting on the reference spacecraft for the range of antenna sizes 
considered. Figures 66 through 68 show the corresponding disturbance 
torques. The curves represent maximum disturbances in LEO and GEO, and 
nominal disturbances in the normal on-orbit orientation in GEO. For a 
given vehicle size these forces differ by a factor of 1-1/2 to 4 and the 
torques by 1-1/2 to 2-1/2 orders of magnitude between LEG and GEO, reflect- 
ing the very much larger disturbances, dominantly due to gravity gradient 
and aerodynamic effects, acting on the vehicle at low orbital altitude. 

The disturbance forces and torques vary by about 2 orders of magnitude 
over the vehicle size range investigated, except for y-axis torques which 
vary by 4 orders of magnitude.* 

Table 11. Scaling Relations for Modular Single Antenna 

(Reference 2) 

Scale Factor = Antenna Diameter = d (meters) 

Range = 15-200 (meters) 

Antenna Power = 1.5 x d (kW) 

Antenna Mass = 0.36 x d^ (kg) 

Array Area Antenna Power x 8.96 (meters^) 

Array Length = 8.85 x Width (meters) 

Array Mass * Antenna Power x 13.41 (kg) 

Mass Total = 135 x d (kg) 

Mass Avionics = Mass Total + Antenna Mass - Array Mass 
Dimensions Avionics Cube =• (Mass Avionics/19)^^^ * L (meters) 

*The very large torques around the y-axis apparently are due to the high 
ratio of antenna to solar panel dimensions in this configuration. Nor- 
mally the torques around the x and z-axes would be expected to dominate. 
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Figure 63. X Axis Forces (From Reference 2) 
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Figure 68. Roll Torques (From Reference 2) 








Intermediate vehicle sizes with antenna diameters of 80 to 100 m 
correspond to the maximum vehicle mass range that can be accommodated by 
a single Shuttle launch, i.e., the range from 12,000 to 18,000 kg (see 
Section 3.4), For these sizes the maximum disturbance torques around the 
X and z axes, in low orbit, are about 20 N-m; those around the y-axis, 
about 100 N-m. 

Based on these data it is apparent that an electric auxiliary thruster 
of any reasonable size, i.e., 0.05 to 0.10 N, cannot provide the maximum 
control torques required in LEO even if the atypical ly large yaw torques 
shown are disregarded. With a long (10 m) moment arm provided by the 
outrigger booms, the thrusters (see Figure 60), used in pairs, would 
not produce control torques greater than typically 1 to 2 N-m. Gimballing 
of the main thrusters would not yield more than about 2 to 5 N-m for 
representative thruster sizes and locations. 

even for much smaller vehicle sizes, electric auxiliary thrusters of 
practical size would not meet control torque requirements in low orbit. 

This control task may therefore require the addition of chemical backup 
thrusters, as previously discussed in Section 3.2.3. 

3.3.3 A uxiliary Thruster Sizing . 

Auxiliary thruster sizing for the selected reference configuration was 
keyed to north-south stationkeeping requirements and the control of geosyn- 
chronous maximum disturbance torques, as discussed in the Boeing report. 
Based on the results obtained in that study, the required thruster size 
ranges from 0.001 N per thruster for a vehicle in the 15-m antenna class 
to 0.01 N per thruster for a 60-m antenna vehicle, assuming a set of six 
thrusters acting together. In the configuration shown in Section 3.2 
(Figure 60), with three thrusters firing in +y or -y direction for N-S 
stationkeeping, the required thruster size would be 0.002 to 0.02 N, and 
still somewhat larger for vehicles in the 80 m antenna class. 

Results of an independent calculation of N-S stationkeeping thrust 
requirements are presented below. 

Let 1/n be the fraction of time each day that is assigned for contin- 
uous N-S stationkeeping thrust operation. (Or equivalently, assume a full 
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day of thrust operation is needed every n'th day.) The thrust acceleration 
a required to maintain such a schedule depends on n and y: 

.-7 


^ = 2.3 X 10 


n 

sin Y 


where y is the arc length on each side of the nodal point over which 
thrusting continues, y = 90 degrees, implies continuous thrusting in north 
and south directions, respectively. The term sin y expresses the loss of 
thrusting efficiency associated with thrust arc extension away from the 
nodes. 

Figure 69 shows a/g versus n with v as a parameter. The limiting case 
is n = 1,Y = 90 degrees, which represents continuous thrusting every day. 
For a 15,000 kg spacecraft the total thrust force reouired under this con- 
dition would be 0.04 N. For n = 2 (thrusting every other day) the *^»"ust 
would have to be increased to 0.08 N. 

Assuming a 3000-second specific impulse and Case-I thruster character- 
istics, the stationkeeping power requirement for continuous thrust opera- 
tion (n = 1) for a 15,000 kg spacecraft mass would be approximately 1.5 KW. 
This compares with 80 tc 100 kW required to operate main thrusters of about 
6 N thrust force, assuming Case-II thrusters at I^p = 1500 seconds. 

3.3.4 Principal On-Orbit aV Requirements 

The principal aV requirements include stationkeeping, repositioning 
maneuvers and the final satellite disposal maneuver to a higher orbital 
altitude. A 10-year life in orbit is assumed. 

3 . 3 . 4 . 1 North-South Stationkeeping 

Data derived from the Boeing study indicate 10-year stationkeeping 
requirements ranging from 552 m/sec for the smallest vehicle of the refer- 
enced category (15 m antenna) and 570 m/sec for an intermediate siz#' vehi- 
cle (60 m antenna), to 613 m/sec for the largest size (200 m antenna). The 
increase is due to radiation pressure compensation requirements. N-S 
stationkeeping alone would require 457 to 518 m/sec depending on the calen- 
dar years of the mission. 

The intermediate requirement above is assumed for propellant mass 
estimates calculated in Section 3.4. Taking a loss in thrust efficiency 
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Figure 69. Thrust Acceleration Required vs Duty Cycle 
In N-S Stationkeeping 




for distributed thrusting about +45 degrees centered at the nodal crossings 
resulted in a N-S aV requirement of 639 m/sec. This requirement corre- 
sponds to those calendar years requiring the most stationkeeping. 

3. 3. 4. 2 East-West Stationkeeping 

For large spacecraft structures with high area-to-mass ratios, A/M, 
such as those considered in this study, the annual E-W stationkeeping aV 
expenditures required to compensate for the erects of earth's triaxiality 
are negligible by comparison. 

Based on data from the Boeing study (Reference 2) and from a NASA 
report by R. R. Lovell and T. A. O'Malley (Reference 12), the E-W station- 
keeping aV required to compensate for solar pressure is shown in Figure 70 
as function of the A/M ratio, with the duty cycle p as parameter. N-S 
stationkeeping requirements are also shown for comparison. The value p * 0 
corresponds to impulsive thrusting, p = 1.0 to continuous thrusting over 
the entire orbit during days when stationkeeping maneuvers are being per- 
formed. In low thrust operation, both E-W and N-S stationkeeping aV 
requirements increase with the length of the duty cycle. The penalty fac- 
tor is ti/ 2 = 1.571 for p = 1.0 relative to impulsive thrust maneuvers 
(where p = 0). Figure 70 shows that the aV required for E-W stationkeeping 

exceeds the N-S stationkeeping aV if the A/M ratio is greater than 
2 

0. 45 m /kg, assuming the same duty cycle is used in either case. 

Typical ranges of A/M values for seven classes of large spacecraft 
investigated in the Boeing study are indicated below in the abscissa in 

Figure 70. The "box structure" class of spacecraft is represented by the 

2 

largest range of A/M values, reaching A/M = 1.67 m /kg for a characteristic 

length L of 1000 m. Such a spacecraft would require a yearly aV of 

180 m/sec for E-W stationkeeping if impulsive thrust maneuvers are assumed, 

1. e., 3.8 times more than for N-S stationkeeping. For the "modular 

antenna" class of spacecraft assumed as typical in the present study, the 

E-W stationkeeping requirements range from 20 to 100 m/sec per year as A/M 

2 

varies between 0.17 and 0.90 m /kg, i.e., for antenna diameters ranging 
from 20 to 140 m. In this case the E-W stationkeeping requirements exceed 
those for N-S stationkeeping if d > 60 m. 
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The above data were shown only in summary form. A more detailed dis- 
cussion, based on equations and parameters that govern stationkeeping 
maneuver requirements, are presented in the following discussion. 

An alysis o f aV Requirements fo r E-W Stat ionke eping 

In Reference 12 the aV requirements for various methods of E-W sta- 
tionkeeping to compensate for solar pressure effects are derived. It is 
shown that circumferential thrust in E or W direction requires only one- 
half the aV expenditure that would be used hy the alternative of radial 
thrusting. It is also shown that this thrust mode requires ISif. more AV 
expenditure than one of direct solar pressure compensation where continuous 
thrust parallel to the sun line would be applied. Considering the nominal 
earth-pointing spacecraft orientation, however, it is apparent that inter- 
mittent thrusting in a fixed E or W direction is much simpler than apply- 
ing thrust components along the three spacecraft axes in periodically 
varying combinations so as to produce a continuously sun-pointing thrust 
vector. The circumferential thrust mode (termed method 2 in Reference 12) 
therefore will be adopted as the nominal mode in the remainder of this 
analysis. 

According to the derivation in Reference 12 the annual aV requirement 
is given by 



(61) 


where 

k = (1 + o) ^ (with 0 average reflectivity of satel 1 ite, '^0.3) 

-fi -1 

S = solar constant, 4.5 x 10 kg m” sec" 

X = mean angular velocity of earth motion around sun, 1.99 x 10"^ 
rad/ sec 

p = duty cycle 

B = eccentricity ratio eVe^ (with e* = max. allowable eccentricity 
and = peak eccentricity that would occur if the initial orbit 
were circular (e^ = 0) and no stationkeeping were applied.) 

The parameter B which defines the desired stationkeeping accuracy can 
also be expressed by the relation 
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= 0.4 


(63) 


whare aL^ is the allowable longitude excursion (in degrees) between sta- 
ticnkeeping maneuvers, (b is a dimensionless quantity.) For very email 
values of aL^, and hence 6, the term (B/sin“^6) in Equation (61) approx- 
imately equals 1.0, and therefore stationkeeping accuracy will have little 
if any influence on aV^ in this case. 

The thruster acceleration required in performing E-W stationkeeping 
is given by 



where 

m = number of days to complete a stationkeeping correction 
0^ = angular velocity of earth's axial rotation, 

= 7.29 x 10”^ rad/sec 

Using results obtained in Reference 12, Figure 71 shows the normalized 

2 

velocity increment AV^/k in (m/sec)/(m /kg) as function of the duty cycle 
p with B as parameter, based on Equation (61), and Figure 72 shows the 
normalized thruster acceleration (m/k) a^ in g's/(m‘'/kg), based on 
Equation (63). 

The curves in Figure 71 show the influence of the principal parameters 

p and B on aV_. However, for representative k values ranging between 0.2 
2 * 

and 2.0 m /kg and for longitude excursions aL^ assumed to be less than 0.1 
degree, 6 is restricted to the range of 0.02 to 0.2, and the aV require- 
ments are chose given by the top curve in Figure 71. For these small B 
values and given values of k, the only way to achieve a reduction in AV 
requirements below the maximum is by reducing the duty cycle p to less 
than 1.0, taking advantage of the fact that aV^ is most sensitive to a p 
reduction in the upper range of its values. 

The multiplicative effect of the A/M ratio, and hence k, on aV^ 
requirements was previously exhibited in Figure 70 for p » 0 and 1.0 assum- 
ing low B values. 
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Table 12 shows results of a parametric evaluation of yearly 
requirements for spacecraft with antenna diameters ranging from 20 to 
140 m. The velocities shown are those for duty cycles p = 1.0, 0.3 and 0 
(impulsive thrust) in circumferential thrusting (method 2) and also the 
corresponding values for direct compensation by thrusting in sun direction 
(method 1). The selected maneuver mode (p = 0.3) requires 34% less aV 
expenditure than operation at p = 1 and 22% less than operation in the 
direct solar pressure compensation mode. The A/M values used in this cal- 
culation are based on data given in Reference 2 for the modular antenna 
class of spacecraft that is considered in this analysis. Accordingly, the 
dependence of A/M on the antenna diameter is given by 

= m {r ^ 

where 

M = 135 X d (kg), (d in meters). \ 
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Figure 71. Normalized Velocity Increment per Year for E-W Station- 
keeping as Function of Duty Cycle with Eccentricity 
Ratio as a Parameter (from Reference 12) 
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Figure 72. Normalized Thruster Acceleration for E-W Station- 
keeping as Function of Duty Cycle with Eccentricity 
Ratio as a Parameter (from Reference 12) 


In addition Table 12 lists nominal spacecraft mass and power and 
assumed stationkeeping thrust, thruster power and acceleration character- 
istics as functions of antenna diameter, and derives values of a^ (m/k) for 
p = 0.3 based on Equation (63) for two stationkeeping accuracy requirements 
aL^ = 0.05 and 0.10 degrees. Finally the number of days (m) is calculated 
during which stationkeeping maneuvers at a 30% duty cycle are to be per- 
formed. Also listed is the stationkeeping cycle duration t^ (in days) for 
the specified aL^ values, t^ is given by the relationship, derived in 
Reference 12, 

t = niLl X 365 ^ X 365 days for e « 1 (65) 

The results show that stationkeeping cycles forAL^ = 0.1 degree range 
from 5.3 to 16.5 days, each cycle requiring 5 days of stationkeeping maneu- 
vers lasting for 7.2 hours per day (30% duty cycle). For aL^ = 0.05 deg, 
the stationkeeping cycles would have to occur twice as often requiring 
7.2 hours of thrusting per day for 2.5 days, but leaving the total annual 
aV expenditure unchanged. 

Mass Variation Effect on Solar Pressure Perturbation 

It should be noted that the aV„ values derived in Table 12 are based 

s 

on fixed values of spacecraft mass. Actually, since solar pressure effects 
depend on A/M and mass varies with propellant expenditure, a more precise 
derivation of aV^ expenditures requires that the variation of A/M be taken 
into account. The iterative technique usee for this purpose is discussed 
in Appendix A. 

3. 3. 4. 3 Repositioning 

A total of five repositioning maneuvers of 180 degrees each is assumed 
as an upper limit in a 10-year mission. For low thrust the maneuver 
requirement is expressed by 

AV = 5.70 ~ (in m/sec) (66) 

^0 

where 

A6 = transfer angle (in degree's) 

t- = time to complete the transfer (in days) 
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Table 12. Characteristics of E-W Stationkeeping Maneuvers 



29 










Allowing 20 days for each transfer, the total aV requirement for five 
repositioning maneuvers becomes 256.5 m/sec. In impulsive maneuvers this 
aV requirement would be reduced by 2. 


3. 3. 4. 4 Orbit Raising Maneuver for Satellite Disposal 


The maneuver requirements are given by 


aV, = V 


aVo = V 



initial impulse 


final impulse 


where 


V 


CO 


r 

r 

r 


0 

1 

2 


= circular velocity at r^, 7908 m/sec 

= Earth radius - 6378 km 
= 6.614 r^ = radius at synchronous altitude 
= radius at altitude of disposed orbit 



(67) 


(68) 


The velocity requirement (aVj + AV 2 ) for disposal to 40,785 km alti- 
tude is 166.5 m/sec. 

3.3.5 Use of Batteries 

The question of whether auxiliary propulsion functions might be sup- 
ported at times (e.g., during eclipse) by the use of batteries was 
investigated. 

Assuming a storage requirement of 10 kW-hr and a voltage of 120 V, 
about one hundred 80-amp-hr cells would be required, at a discharge volt- 
age of 1.25 V/cell. Depth-of-charge of 8 O 1 ; is acceptable for this appli- 
cation since only 60 charge/discharge cycles per year are required. Total 
losses of 20 % was assumed, 10 % each during charge and discharge operation. 
The battery set would have an estimated mass of 450 kg including battery 
charger. 
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The use of batteries for this purpose does not seem weight efficient 
since chemical backup thrusters will be available in any case, as discussed 
previously, and can substitute for electric auxiliary thrusters when power 
is unavailable. Secondly, momentum storage devices, being carried to 
handle routine ACS operations, would obviate the use of electric ACS 
thrusters during the 72-minute eclipse periods occurring at the spring and 
fall equinox seasons. 

The use of batteries as a buffer against peak power demands by the 
eclipse APS, for example to avoid a conflict with payload power require- 
ments, does not appear necessary. Power sharing in the orbital mission 
phase often involves under-utilization of the available power, rather than 
a need for managed power allocation under peak loads as is discussed in 
Section 3.4. In missions with unusually large payload power requirements 
this question may require further study. 

3.3.6 Thruster Cant Angle Requirement 

Auxiliary thruster location and orientation generally require detailed 
design trades depending on overall spacecraft configuration and mission 
profiles. Significant constraints will be imposed on thruster placement 
and orientation to avoid thruster plume impingement on the arrays, espe- 
cially for spacecraft with large solar arrays, as discussed in Section 3.2. 

Referring to Figure 60, the boom-mounted thrusters may require an 
offset or cant angle large enough to avoid plume impingement on the nearest 
solar panel areas under worst-case panel orientation which will occur twice 
each day. To define the best arrangement a trade is required between the 
propellant pe".*lty due to the cosine-loss resulting from canting the 
thrusters versus the mass of added boom structure that would be required to 
reduce the cant angle. 

3.3.7 Gravity Gradient Strategy 

During the ascent phase of the mission, gravity gradient torques can 
be employed to assist In orienting the spacecraft for optimum array Illumi- 
nation. The spacecraft attitude maneuvers required are summarized In 
Table 13. In the atmospheric drag region, the solar arrays are feathered 
for minimum drag. The spacecraft is either maintained fixed at ♦ (roll 
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Table 13. Spacecraft Attitude Maneuvers for Optimum Sun Alignment 
During Primary Thrust Phases 


1, Lowest Altitude 

Roll Only 

Excursions ±38* for Sun Angle * 52* 

(ARRAY FEATHERED FOR 


t75* for Sun Angle » 15* 

MINIMUM DRAG) 


Fixed 90* for Sun Angle < 15* 

2 , Low- to- 1 ntermed 1 a te 

Roll Only 

SAME AS ABOVE 

Altitudes 



(ARRAY ARTICULATED) 



3. Intermedlate-to-HIgh 

Roll 

Modified Periodic Excursion In Same 

Altitudes 


Range 

(ARRAY ARTICULATED; 

Yaw 

Excursions to t45* GENERALLY NOT IN 

PLANE CHANGE SIMUL- 


PHASE WITH ROLL EXCURSIONS. DEPEND- 

TANEOUSLY WITH ORBIT 


ING ON NODAL POSITION RELAT.IVE TO 

RAISING) 


SUN LINE 


angle) = 90 deg or rolled from ±38 to ±75 degrees on sun angle (see 
Section 2.6.1 for a discussion of seasonal sun angle variation rela+’ive to 
the orbit plane). Gravity gradient is useful in providing restoring 
torques towards 4i = 0. The maximum gravity gradient torques are near 
(j> = ±45 deg, which are consistent with maneuver requirements. The best 
utilization of gravity gradient occurs when the period of roll motion is 
close to that of the orbital period. At low to intermediate altitudes, 
above the drag region, the solar arrays are rotated about their axis, and 
the roll motion is the same as described above. At intermediate to high 
altitudes, when combined orbit boosting and plane change are implemented, 
the roll motion is similar to that described above, but is modified to 
accommodate additional yaw motion for optimum solar array Illumination. 

The gravity gradient assist is most useful at low altitudes, where 
torque requirements are greatest, for reducing the size of inertial momen> 
turn storage devices in the attitude control system. Since the roll motion 
is driven primarily by gravity gradient, rather than by spacecraft attitude 
control moment gyroscopes or reaction control, solar array deformations are 
also reduced. If possible, spacecraft mass distribution should be chosen 
to permit gravity gradient operation slightly below the natural oscillation 
frequency. Thus, inertial design of the spacecraft would be constrained, 
but large roll angles could be maintained with small driving torques. 
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3.4 INTEGRATED PROPULSION ADVANTAGES 


The advantages of integrated electric propulsion are assessed by 
weight and concommitarit transportation cost comparisons with chemical 
propulsion. On-orbit power sharing with electric propulsion is another 
advantage examined for the reference generic spacecraft. 

3.4.1 Comparison of Electric and Chemical APS 

The aV requirements for APS were presented in Section 3.3.4. They 
are summarized here in Table 14, which lists requirements other than those 
for E-W stationkeeping. The iterative procedure described in Appendix A 
was used for determining aV^, the E-W stationkeeping requirements. 

Comparison was made between a chemical APS operating at 300 seconds 
I^p, and an electric APS at 3000 seconds I^p. As discussed in Sec- 
tion 3.1.1, repositioning and disposal in the electric case were performed 
by the main thrusters (in the integrated EPS) at 1500 seconds I^p. 

Table 15 lists aV requirements and resulting mass characteristics for the 
electric and chemical APS cases for spacecraft with 20, 60, and 100 m 
antenna diameters as well as the diffe>"erces in total mass launched. The 
mass characteristics are shown graphically in Figure 73 for comparison. 

Note, that for the 60 m antenna spacecraft, the use of electric APS 

★ 

saves over 7000 kg of total launched weight. At the Shuttle cargo limit 
of 29,500 kg, electric APS affords about a 50% increase in reference 
spacecraft mass. 

In Table 15 the total mass values and + FSE are derived from 
the initial mass in orbit Myj by adding the propellant mass consumed 
during low thrust ascent from LEO to GEO (aV » 6015 m/sec) and the mass ot 
mission-peculiar flight support equipment (FSE), such as payload cradles, 
which is assumed to be 5% of the gross mass Myj^. The primary electric pro- 
pulsion system mass including tankage is included in the reference mass 
The auxiliary propulsion system dry mass is included in 


1 * 

The terms weight and mass are used interchangeably in this discussion. 
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Table 14. Auxiliary Propulsion aV Requirements 
Other Than for E-W Stationkeeping 
(10-Year Mission), in m/sec 


ELECTRIC APS 

1 . N-S Stationkeeping 

(Assuming + 45® thrusting arc) 

- Low Value 

564 

- High Value 

639* 

2. Repositioning (5 Events) 

256.5 

3. Disposal to 40,785 km Altitude 

■;66.5 

CHEMICAL APS 

1 . N-S Stationkeeping 

- Low Value 

508 

- High Value 

575* 

2. Repositioning (5 Events) 

128.3 

3. Disposal to 40,785 km Altitude 

166.5 


*Th1s value used in subsequent analysis 

3.4.2 Launch Cost Savings Achievable with Inte grated Electric Auxiliary 
Propulsion 

NASA's Shuttle Launch cost reimbursement policy (Reference 13) defines 
weight dependent and length dependent user charges, respectively, as 

<^W “ '-“S X " (<teHcated launch cost) 

'^L ' ‘ ”5 « to? l frl F 9 i°caiac ;^ * launch cost) 
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Table 15. &V and Mass Characteristics of Electric and Chemical APS Spacecraft 
(Reference Configuration: Modular Antenna System) 
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where dedicated launch cost is given as $ 30.2 million (1981 dollars). 

The user will be charged the larger of the two cost figures. The full 
dedicated launch cost will be charged if either spacecraft weight or 
length exceed 75% of full capacity, i.e., 22,110 kg or 13.72 m, 
respectively. 

Figure 74 illustrates this reimbursement policy in terms of cost con- 
tours plotted in a mass versus length diagram. The graph shows the fixed 
cost plateau reached when jiiass or length exceed the 75% of capacity level. 
The dashed diagonal designates the breakeven points of weight and length 
dependent charges. Using the dedicated Shuttle launch cost, the weight 
and length dependent charges are 

Cy = $1,384 million per 1000 kg 

C|^ = $2.30 million per meter 
The slope of the breakeven line is defined by 

M = 2.30 = 1.611 X 10^ kg/m 

Lg^ 1.384 

If lauru,h charges were always weight-dependent and varying linearly, 
the launch cost savings due to weight reduction would be defined simply by 
the relation 

Acost -- $1,384 aM (millions per 1000 kg) 

For 3xample, the 7,209 kg weight savings for the nominal 8,100 kg 
spacecraft mass (see Table 15) would translate into a launch cost saving 
of approximately $10 million. 

Actually, under certain conditions these savings will not always be 
fully realized under the NASA launch cost reimbursement policy, for exam- 
ple, if spacecraft dimensions make the launch cost length-dependent. 

Figure 75 shows bars of weight savings aH for several values of assumed 
spacecraft length. In one case, the entire aM is in the weight-critical 
region, in the second case it straddles the cost breakeven line aru in 
the third case it is entirely in the length critical region of the cost 
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Figure 75. Examples of Launch Cost Savings 
Variation wUh Cargo Length 
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contour diagram. Accordingly, the cost savings reflect either all or part 
of the weight savings; or, in the third case, no cost savings are realized 
at all in spite of large weight savings. 

Another effect of NASA's nonproportional cost allocation is the partial 
loss of potential savings if the gross weight exceeds the 75% limit. 

Significant launch cost savings with integrated electric propulsion 
are anticipated for weight dependent payload configurations. Length depen- 
deiit configurations will also benefit because of the reduced volume required 
for propellant storage when compared with chemical propulsion. 

3.4.3 Power Sharing in Orbit 

Shared use of available solar array power by the integrated electric 
APS and the spacecraft housekeeping systems and payload(s) during the 
orbital mission phase tends to increase overall system cost-effectiveness, 
a further incentive for adopting the integrated APS concept. 

To characterize the degree of power utilization during the mission, 
the ratio of power used on-orbit to power used in ascent is introduced. 

This ratio, r, serves as a measure of power utilization effectiveness. 

Figure 76 shows contours of constant r in a graph of on-orbit versus ascent 
power requirements on a logarithmic scale. 

Consider as an example a mission which requires 200 kW initially for 
orbit-transfer by primary electric propulsion. After a power loss of about 
40% during the transfer phase due to solar array radiation damage, the 
remaining power level .wuld be 120 kW. On orbit power required by the pay- 
load and housekeeping subsystems is assumed to be 30 kW, i e., calling fcr 
only 25% utilization of available power (Case 1). 

By using auxiliary electric propulsion with a nominal power allocation 
of 30 kW for on-orbit maneuvers the utilization effectiveness is raised to 
50% (Case 2A). At times of primary propulsion use for major on orbit 
maneuvers, however, all available power is allocated to the EPS, thus 
raising the utilization effectiveness to 100% (Case 2B). 

Considering the launch weight reduction achievable through use of the 
integrated electrical APS as discussed in Section 3.4.1, the solar array 
size could actually be reduced by 30 to 35% without reducing the initial 
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Power Utilization Examples 





thrust acceleration. The solar array power would be 130 kW initially, the 
degraded power level about 80 kW. Assuming the same payload and auxiliary 
propulsion power requirements on orbit as before, the utilization factor 
becomes x 100 = 75% (Case 3A) with an occasional increase to 100 per- 
cent during major maneuvers (Case 38). 

In summary, the integrated APS concept will increase power utiliza- 
tion by a factor of 3 (or 4) in Case 3, while also permitting a one-third 
reduction in launch weight, solar array and primary electric propulsion 
size. 

In Case 2 where the spacecraft uses the same solar array and primary 
electric propulsion system as in Case 1 with a one-third reduction in 
launch weight, the time required for orbital ascent would be substantially 
less than in Cases 1 or 3. The power utilization factor would be raised 
by 2 (or 4) compared with Case 
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4. DISCUSSION OF RESULTS 


This study has been primarily concerned with the LEO to GEO mission. 
Several results from the study, such as solar cell cover thickness optimi- 
zation, optimum I^p range, importance of thruster efficiency, and low 
altitude flight apply to orbit transfer portions of the analysis. They 
are discussed separately below. Next, a specific mission has been selec- 
ted to illustrate how study results can be applied to particular configura- 
tions of interest. Results from both orbit transfer and integrated pro- 
pulsion analysis are used to define the integrated EPS. 

4.1 ORBIT TRANSFER 

4.1.1 Solar Cell Cover Thickness 

It was found that, for the LEO to GEO transfer mission, there appea>^s 
to be no compelling reason to reduce the thickness of solar cell covei^s 
much below 6 mils. The calculations were made for an assumed payload 
fraction of 33%. For the lowest assumed thruster efficiency function the 
cover thickness for minimum t"*ansfer time ranged from 5.5 mils to 6.5 mils 
as the specific mass a varied from 0.01 to 0.02 (W/kg). For the highest 
assumed thruster efficiency function the corresponding cover thicknesses 
ranged from 3-1/2 to 4-1/2 mils for specific masses varying from 0.005 to 
0.01 (W/kg). For the shortest transfer time case (highest efficiency and 
lowest a), the transfer time for 6-mil covers was only 12% larger than the 
minimum time for the 3-1/2 mil covers. However, the corresponding 
residual power on orbit was 15% greater for the 6-mil covers. 

Tradeoff calculations for a given mission would be needed to define 
an optimum thickness for a specific application. However the weak 
dependence of orbit transfer time on cover thickness, i.e., the f’atness 
of the minimum, shows that cover thickness is not a sensitive parameter so 
that 6 mils is a near-optimum thickness for a wide range of parameters. 
Improved solar array performance might be expected from thinner covers for 
missions below the Van Allen belt; but for the LEO to GEO application, 

6-mil covers appear to be adequate. 
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4.1.2 Optimum Up Range 

The optimum I^p for transfer from LEO to GEO is the value that gives 
the shortest transfer time for a given efficiency. It depends on payload 
fraction but is essentially independent of both the specific mass param- 
eter a and the efficiency. Its value ranges from 2600 seconds at 50* 
payload down to 1400 seconds at 20* payload. 

The functions relating n to I^p for a fixed transfer time are fairly 
flat near the optimum I^p. Because n generally increases with I^p for 
electric thrusters, the I^p for minimum transfer time for a given type 
thruster occurs at a higher value than the optimum I^p, the difference 
being larger for the less efficient thrusters than for the more efficient 
ones. For example at 33* payload, the least efficient thruster had an 
Ijp of 4000 seconds for minimum transfer time as compared to 2100 seconds 
for the most efficient thruster. 

In the neighborhood of the minimum transfer time (when the slope of 
the efficiency function is equal to that of a constant transfer time line), 
variation of I^p has very little effect on transfer time. But it does have 
a strong effect on the trade between dry propulsion mass and propellant 
mass. Reducing the I^p only slightly results in a large reduction in the 
more expensive dry propulsion mass that is compensated by a nearly equal 
gain in the less expensive propellant mass. For example, using the inter- 
mediate efficiency relation at 33* payload, the minimum transfer time is 
174 days for an I^p of 2300 seconds. The corresponding values of dry 
propulsion mass and propellant mass are 1.06 and 0.80 times the payload 
mass. If the I^p were reduced to the optimum value of 1750 seconds, for 
this case, these values would become 0.84 for dry propulsion and 1.14 for 
propellant. The transfer time, however, would be increased from 174 days 
to 190 days, or about 9*. 

A reasonable compromise is to select an I^p midway between the opti- 
mum value and the minimum transfer time value, about 2000 seconds for the 
present example. The transfer time is then about 180 days, the dry pro- 
pulsion mass is 0.96 and the propellant mass is 0.92 times the payload mass. 
Using this appraoch, the I^p ranges from 2850 seconds at 50* payload to 
1700 seconds at 20* payload for the intermediate efficiency relation. 
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Higher values of I^p are encountered for planetary missions that 
require much greater total impulse than that needed for LEO to GEO trans- 
fer. Lower values are optimum for lower impulse missions such as small 
changes in orbit altitude ranging from the Shuttle parking orbit to the 
beginning of the radiation belt for spacecraft using integrat®ci propulsion. 
For these missions high energy chemical rockets and resistojets become 
attractive. However, for the total impulse range of LEO to GEO missions, 
the useful I^p range is roughly from 1500 to 2500 seconds. This may also 
be true for multimission spacecraft, such as an electric tug, making 
several low altitude orbit changes. This aspect of Shuttle capability 
extension needs further study. 

4.1.3 Importance of Efficiency and Specific Mass 

The two most important thruster system parameters affecting transfer 
time from LEO to GEO are the thruster efficiency (the ratio of ideal beam 
kinetic power to input electric power) and the thruster specific mass 
(mass per unit input electric power of the dry propulsion system). The 
importance of efficiency is illustrated by the 33% payload example having 
a specific mass of 0.024 kg/W. For the lowest efficiency function 
assumed for the study the minimum transfer time is 272 days at 4000 seconds 
Ijp and 48% efficiency. For the intermediate efficiency function the 
time is reduced to 174 days at 2300 seconds I^p and 51% efficiency, ^or 
the most efficient thruster the time is only 147 days at 2100 seconds I^p 
and 59% efficiency. These values are shown in Table 16. 

The Case I and Case II efficiencies are based on present data; the 
Case III function is an estimated performance that may be attainable in 
the near future. It is not so much the efficiency differences that are 
important as it is the combination of increased efficiency with reduced 
Ijp. For instance, the Case II example is only 3% higher efficiency than 
Case I, yet the transfer time is reduced by 98 days as a consequence of 
the Up reduction from 4000 to 2300 seconds. 

A slightly different example. Case II efficiency at optimum Ijp and 
33% payload, is used to show the Importance of the specific mass parameter 
o on transfer time. Table 17 compares results for o ranging from 0.016 to 
0.032 (kg/W) for this case. 



Table 16. Effect of Thruster Efficiency on Transfer Time 
(33% Payload, a = 0.024 kg/W) 




Table 17. Effect of m on Transfer Time 

(33% Payload, n = 45% at 1750 seconds I^p) 


a 

Transfer Time 

(kg/W) 

(Days) 

0.016 

120 

0.024 

193 

0.032 

259 


! 

• « 


4 

*, 

P 

i 


The specific mass parameter is clearly of great importance. The 
transfer time is reduced 66 days by reducing a from 0.032 to 0.024 (kg/W). 

For the LEO to GEO transfer mission, thruster technology effort should 
be focused on increasing efficiency at moderate I^p (1500 to 2500 seconds) 
and reducing the specific mass as much as possible. Any thruster capable 
of achieving both of these objectives simultaneously would be very attrac- 
tive for this .. ■ jSion. 

4.1.4 Fast Orbit Transfer 

In addition to reduction of LEO to GEO transfer time through efficiency 
improvement and specific mass reduction, it is always possible to achieve 
further reductions by giving up payload. At some point this defeats the 
primary advantage of electric propulsion which is its ability to transport 
larger payloads than can be handled by chemical propulsion. Table 18 
summarizes the trade between payload and transfer time for the low speci- 
fic mass of 0.012 kg/W and the highest efficiency function assumed for an 
electric thruster. Transfer times for tfe resistojet, operating at 
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Table 18. Effect of Payload on Transfer Time for Electric and 
Resistojet Thrusters 


Payload Fraction 

W 

Transfer Time (Days) | 

Electric 

Resistojet 

10 

31 

18 

20 

39 

32 

30 

52 

118 

40 

75 

- 

50 

123 

- 


90% efficiency and 900 seconds I^p are also shown for comparison of this 
thruster type. 

This table further illustrates the importance of high efficiency at 
moderate ^P’ combined with low specific mass. It also show that the resisto 
jet may not be suitable for the integrated LEO to GEO mission. It shows an 
advantage only for small payload fractions, for the assumed I^p of 900 sec- 
onds. It would quickly lose this advantage if the I^p were reduced even 
slightly below this value. For example, at 800 seconds I^p and 90% effi- 
ciency, the time would be increased from 32 days to 60 days for 20% 
payload. 

4.1.5 Low Altitude Flight 

It was found during the course of this study that solar electric 
propulsion could be extended to altitudes well below the optimum Shuttle 
parking orbit. For an assumed specific mass of 0.024 (kg/W), n “ 42% at 
1500 seconds I^p and an area density of 30 kg/m a spacecraft with 60% 
payload could maintain altitude at 150 km. This is done by feathering the 
solar array for minimum drag (oriented edge-on to the flight direction) 
and rolling it for maximum illumination to produce an average occulted 
power of about 40% of the full array power. The greatly increased ratio 
of solar array power to drag area opens up the possibility of continuous 
flight in the aerodynamic drag range of altitudes that extend from 500 km 
down to 150 km. Previously the lower end of this range was considered to 
be unacceptable for solar electric propulsion (Reference 14). 
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A comparison of the feathered flight mode with a nuclear power source 
shows comparable performance down to 150 km. 

4.2 SPECIFIC MISSION EXAMPlE 

The specific mission choset. to illustrate application of integrated 
electric propulsion systems is summarized in Table 19. It is a modular 
antenna spacecraft with a 60 m diameter antenna. Its reference configura- 
tion is illustrated in Figure 60 (page 102). The mission is to operate the 
spacecraft for 10 years in geosynchronous equatorial orbit after having trans 
Terred it from Shuttle's parking orbit at 250 km altitude and 28.5 degrees 
inclination. Stationkeeping requirements on-orbit are ±0.1 degree. Upon 
conclusion of on-orbit operations, the spacecraft is transferred to higher 
altitude for disposal. 

The APS AV requirements for this mission are sunwiarized in Table 20. 

N-S stationkeeping requirements are derived from Section 3.3.4. 1. E-W 
stationkeeping requirements are taken from Table 12 in Section 3. 3. 4. 2 
(page 129) for Method 2, o = 0.3. The E-W requirements are several percent 
high because the effect of spacecraft mass variation, discussed in Appendix 
A, was not taken into account. Disposal requirements are given in Section 
3. 3.4.4 (page 130). A contingency of 20% has been added to account for 
momentum wheel unloading during attitude control operations, and for pro- 
pellant residuals. Since the total AV on-orbit is very close to the value 
used in Section 2.8.1 (see Table 5, page 57), the results shown in that 
analysis may be used directly. 

The beginning-of-life power can be estimated from the anticipated re- 
sidual on-orbit power shown in Figure 42 (page 76). A residual power in 
excess of 50% is anticipated for the specific case chosen. The residual 
power requirement of 90 kW (Table 19) thus results in a value of 180 kW BOL. 
Using a specific mass, * 0.024 kg/W, the CPS dry mass is (0.024)(180) 
(1000) » 4320 kg. The payload mass is then the spacecraft dry mass from 
Table 19 less the EPS dry mass. Hence, the payload mass * 13,500 - 4320 » 
9180 kg. 

For a thruster exhibiting Case II efficiency characteristics, the 
tradeoff curves in Figures 27 through 30 (pages 59-62) apply to this ex- 
ample. At 33% payload fraction, for instance. Figure 28 is applicable. 
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Table 19. Specific Mission Requirements 


Configuration 

Modular Antenna 
(see Figure 60) 

Antenna Diameter 

60 m 

Spacecraft Dry Mass 

13,500 kg 

Initial Altitude 

250 km 

Final Operation Altitude 

35786 km 

Initial Inclination 

28.5“ 

Final Operation Inclination 

0° (equatorial orbit) 

Disposal Altitude 

40786 km 

Lifetime On-Orbit 

10 years 

On-Orbit Stationkeeping Accuracy 

±0.1“ 

On-Orbit Power 

90 kw 


Table 20. aV Requirements for APS, 60 m Modular 
Antenna Spacecraft 


Maneuver 

aV (m/sec) 

North-South Stationkeeping 

639 

East-West Stationkeepinq 

545 

Disposal 

167 

Sub-total 

1451 

20% Contingency 

290 

Total 

1741 


Minimum transfer time occurs at 2300 I^p (see Table 6, page 64), but it is 
better to operate at slightly lower I^p for the reasons given in Section 
2.8.1. Selecting 2000 sec as the operating point yields a thruster effi-. 
ciency of 47. 5r, and 187 days transfer time at an electric propulsion system 
mass * (2.0) (9180) = 18,360 kg. The total spacecraft mass launched in the 
Shuttle cargo bay (excluding flight support equipment) is then 18,360 + 

9180 » 27,540 kg. The BOL thrust level for the primary thrusters is com- 
puted from Equation (7) to be 8.7 N. 
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The APS thrusters operate at 3000 sec I^p and approximately 0.1 N 
thrust as discussed in Section 3.3.3. At Case I efficiency, the APS 
thrusters consume 4.0 kW per Equation (7). Checking on-orbit power capabil- 
ity, Figure 42 (page 76) shows 53% residual power at 2000 sec I^p close ■ 33% 
payload ratio. This corresponds to 95.4 kW at beginning of orbital opera- 
tions, which is sufficient to support the payload plus the APS thrusters. 

The EPS characteristics for the 60 m mod’jiar antenna spacecraft are 
summarized in Table 21. 

Table 21. EPS Characteristics for 60 m Modular Antenna Spacecraft 


BOL Power 

180 kW 

Mass Properties 

Dry Weight 

4320 kg 

Propellant Weight 

14040 kg 

Total Weight 

18360 kg 

Orbit Boosting Thrusters 

BOL Thrust 

8.7 N 

Specific Impulse 

2000 sec 

Efficiency 

47.5% 

APS Stationkeeping Thrusters 

Thrust 

0.1 N 

Specific Impulse 

3000 sec 

Efficiency 

37% 

Power 

4.0 kW 

LEO to GEO transfer time 

187 days 

Orbital lifetime 

10 years 
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5. CONCLUSIONS 


(1) For a wide range of ei raft and nission parameters, 6-mil 

covers are near optiipum ror LEO to GEO transfer. 

(2) The most favorable I$p range for primary propulsion for this mis- 
sion is roughly from 1500 to 2500 seconds. 

(3) High efficiency in this I$p range, and low specific mass are 
very important for short transfer times. 

(4) Fairly short transfer times are feasible with near term tech- 
nology. With demonstrated thruster efficiency, the transfer 
time is about 175 days for 33% payload. 

(5) Very short transfer times may be possible with further tech- 
nology advances. For the assumed Case II thruster efficiency 
and a = 0.012 kg/W, the transfer time for 30% payload is 52 days. 

(6) The resistojet is only attractive at low payload fractions for 
this mission. 

(7) With the feathered solar array, rolled for optimum illumination, 
continuous flight is possible at altitudes above 150 km. 

(8) Substantial residual power is available ^or on-orbit power 
sharing with payload. At 33% payload, 57% of BOL power is 
available. 

(9) Large propellant mass savings, through integrated electric APS, 
results in 50% increase of maximum Shuttle launched reference 
spacecraft mass. These savings ar'e realized principally through 
reduction in auxiliary propellant mass and the resoHing, e'en 
larger, reduction in primary propellant mass requirements. 

(10) The launch weight reduction from integrated electric APS can lead 
to Shuttle launch cost savings as large as $10 million for a 
spacecraft of 8100 kg net mass studied herein. However, the 
results are strongly influenced by the spacecraft stowed length 
under NASA's Shuttle launch cost reimbursement policy. 

(11) Electric APS and primary EPS thrust vectoring meet attitude con- 
trol requirements. Chemical auxiliary thrusters only provide 

a backup function. 

(12) Gravity gradient torques can be used to sustain spacecraft roll 
maneuvers for optimum solar panel illumination during primary 
thrust phases of the mission. Operation near resonance should 
greatly reduce control moment gyroscope requirements. 
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(13) With the initial gross weight reduction achieved through 
integrated electric APS, the solar array and primary EPS used 
for LEO to GEO transfer Cun be significantly reduced in size 
while still producing a similar acceleration as before. 
Alternatively, if not reduced in size, the system would achieve 
a faster transfer due to increased acieleration. 

(14) East-west stationkeeping aV due to solar pressure on large 
spacecraft strucvures can be three to four times greater than 
north-south staticnkeeping aV. There is practically no effect 
of east-west staticnkeeping accuracy on aV requirements. 



APPENDIX A 


EAST-WEST STATIONKEEPING AV REQUIREMENTS INCLUDING THE EFFECT OF 
SPACECRAFT MASS VARIATION (ITERATIVE PROCEDURE) 


Spacecraft mass variation due to stationkeeping propellant expenditures 
are relatively minor, i.e., only a few percent, in the case of electric 
auxiliary propulsion. However for an accurate performance comparison with 
chemical auxiliary propulsion maneuvers where such propellant expenditures 
will be very much larger, the effect of mass variation of solar pressure 
perturbation must be taken into account. 

Using the aV^ requirements derived in Section 3. 3. 4. 2 (Table 12) as a 
first step in an iterative process, a good approximation of actual aV^ 
requirements and the corresponding propellant expenditures is obtained 
by the algorithm derived below. 


We use the average spacecraft mass defined by 


M = M„ + 4 Mod = 
ave 0 2 PP o 


1 . I 




in approximating an average value of k 


(A-1) 


k.„^ = 1 .3 rr 

ave M 


ave 


in Equation (61) which derives aV^ as a function of k or A/M. 

The term aV used in Equation (A-1) is the total aV expended on orbit, 
for stationkeeping and other maneuvers, av^ is only a part of this total 
aV expenditure and must be separated from the other part designated by 
In the equa.ion relating AV^ to 


AV^ = Q 


•^ave 


M. 


ave 


(A-2) 
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the terms aV^ and aV^ used in defining must be segregated as follows 


aVs = 1.3 


n 

M 


(1 




+ aV 




(A-3) 


The coefficient Q used as a simplifying notation in Equation (A-2) and 
(A-3) is the product of all factors except k = A/M that appear on the right 
side of Equation (61). 


To further simplify Equation (A-3) v/e substitute the terms 


e^Vr/Ispg = 

e^Vs/IspQ = 5^ 

l-3Q^=iV,0 

where the quantity aV^, is the aV^ value that was previously obtained by 
assuming a constant mass M^ (Table 12). With these substitutions Equation 
(A-3) is reduced to 

(1 + S^) = 2 aV^^ (A-4) 

or 

aV^ (1 + = 2 aV^q (A-5) 


From this equation the corrected values of aV^ can be readily derived given 
V ^''so- 

Note that the principal approximations made in this procedure are the 
assumptions that the average mass is the mean between the masses at 

the beginning and end of the orbital mission phase and that is pro- 
portional to 1/M . The approximations are better for electric propulsion 

where the propellant expenditure is relatively minor than for the chemical 
APS case. However, for purposes of this analysis the approximations inher- 
ent in the above approach are considered adequate. 

Note also that a differentiation of I^p values between those maneuvers 
performed by the main and the auxiliary electric thrusters is approximate 
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and should be made in calculating the corrected values using an equa- 
tion similar to (A-5). In these calculations we assumed specific impulse 
values of 1500 seconds for the main and 3000 seconds for the auxiliary 
electric thrusters vs 300 seconds for chemical thrusters in auxiliary pro- 
pulsion system performance comparisons. Table 14 lists AV requirements 
other than those for E-W stationkeeping, i.e., the partial tM requirements 
adding up to aV^ in Equation (A-3) to be used as input terms in solving 
for aV^. 

Table 22 gives the results of relating corrected and initial aV^ 
values for electric and chemical propulsion maneuvers, based on the above 
I^p data. (In solving the transcendental velocity equation (A-5) for the 
unknown variable, aV^ is assumed as input and AV^^ obtained as output 
quantity.) The results are plotted in Figure 77 which shows that in the 
range of interest the corrected aV^ values for electric propulsion are 


Table 22 . Solution of Equation A-5 for Iterative Correction of aV 

(Velocities in m/sec) 


AV^ 

Ss 

^''ns 

^NS 

AV^ 


so 

ELECTRIC APS 
(I 5 P = 3000 s) 

(3000 s) 

(1500 s) 


300 

1.0102 

639 

1.022 

423 

1.02 ^ 

309 

600 

1.0206 

639 

1.022 

423 

1 . 02 a. 

621 

900 

1.0311 

639 

1.022 

423 

1.0292 

937 

1300 

1.0452 

639 

1.022 

423 

1 .0292 

1363 

— 

CHEMICAL APS (I^p * 

300 sec) 





300 

1.1074 

870 

1.344 



367 

600 

1.2264 

870 

1.344 



794 

900 

1.3581 

870 

1.344 



1271 

1300 

1.5561 

870 

1.344 



2009 


155 


H 


f 




ORIGINAL PACE IS 
OF POOR QUALITY 



Figure 77, Iteratively Corrected Values aV Versus aV (10 Year Mission) 



3 to 5 percent less than the initial values while for chemical pro- 

pulsion the differences can be as large as 35%. As previously explained, 
these reductions are due to the use of average A/M ratios in the cal- 
culation rather than the initially assumed higher A/M^ ratios that did not 
take propellant mass variations into account. Clearly, this effect is 
much larger in the chemical propulsion case. 
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